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r.HAPTER-I 


A smrnY of SOME STOCHASTIC MODELS ON FERTILITY ANALYSI S 


INTRODUCTION 


1.1 Fertility analysis is a very important item of 
investigation in population studies. It refers to the 

additive component of the growth of population. 

Stochastic modelling, of late, has been proved to 

be a very useful technique for fertility investigation. 

This is because several socio-economic and biological 

variables affecU the fertility pattern are either 

unobservables or their measurements are imprecise due 

to non-sampling errors. For example, the estimation of 
proportion of biologically sterile females in a population 
is beset with many problems baffling the exact estimation 
of its size. proportion of mothers without children cannot 

; :be considered as the relevant data for the estimation because 
there are females who remain childless without being sterile 
I under the constraints of physical, social, personal and 

||iychologlcal factors. Under the situation, a judicious 
llse of statistical models iox the description of the process 
I’ahd’ estimation of the relevant parameters often has been 

’ II ' 'I. 

; |l found to bo indlspousoblo^ ■ , - • , 

■ ^^%->:''iFertlllt^r, analysis stochastic models. ^ ^ of. late. ^ 

'-'k* .'k v ■ ’ ' 




Demographers; especially because of certain special aavautcLg^*. 
of the same over the traditional deterministic and actuarial 
approaches. 0^ of the major advantages of stochastic 
models in fertility analysis is to derive small sample estimates 
regarding various fertility parameters with exact specification 
of sampling errors. Also to evaluate the status of various 
possible estimators of the Demographic parameters for their 
relative acceptability under the given set up of a population. 
In fact, the same not generally available necessary in 
actuarial approaches as the latter has the advantage of 
a large number of case records. Besides this, certain 
demographic parameters are basically unobservables but 

playing very important role while characterizing qualitatively 
and quantitatively the structure of a population. Let 

us take the : example of estimating the parameter concerning 
the proportion . of Biologically sterile couples in a population. 
It may b^^ noted that Biologically sterile couples are 

unidentifiabil' from the traditional data of sample surveys. 
Because the|J^al way to idauUty them viz. couples who 
did not ha^W children during the whole reproductive 

no means, provide any clue relating 
sterile couples. A sizeable 


(or Marital) -v span by 

* "l' 

to the size^fofv the Biologically 
number of ^'wuples do appear In records who did not have 
any survii|n|:; chiidr^ biologically not sterile. 

i^Were "are...a^.^number of couples (especially 

Western Countries) who intend 









to remain childless although not being biologically sterile. 

The estimation of the proportion of biologically sterile 

couples, is thus not feasible without appropriate stochastic 
modelling. Second example relating to the effecacies of 
stochastic modelling in fertility analysis arises in the 

case of estimation of Post Partum Amenorrhoea. Estimation 
of post partum amenorrhoea (P.P.A) is beset with many 
practical difficulties (Biswas (1963), Biswas (1967)1- 

Two different types of non-sampling errors do vitiate the 

estimation of P.P.A. and these two errors operate in reverse 

directions in under-estimating and overestimating respectively 
the P.P.A. Besides that, there exists heterogeneity in 

the data of P.P.A. v;hich makes the estimation of P.P.A. 

still difficult -, as the modelling of the P.P.A reveals 

that the same conforms to a mixture of two independent 

random variables (Yadava (1966), Biswas (1971), Singh 
and Bhaduri (1971)). A practical way out to this difficulty 

using stochastic modelling approach is to' assume certain 
plausible probability distribution for describing the fecundable 
period, gestation period and the P.P.A. (Talwar (1966)). 
If the distributions are independent the convolutions of 

the three would be the distribution between two consecutive 

order of conceptions (or births) on the assumption of one 
to one correspondence between the conception and the birth. 
Hence the synthesized (convoluted) model may be fitted 
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in the data corresponding to the waiting time between 
two consecutive births (or conceptions)V The reporting 
of which is much less errorneous than that of the distribution 


of 

the post 

partum 

amenorrhoea. 

The 

fitted data 

in 

the 

model 

provide 

estimates of 

the 

parameters in 

the 


synthesized distribution of which one or two may be pertinent 
to the period of amenorrhoea- The example illustrates the 
applicability of stochastic modelling of estimating the P.P.A. 
without drawing any data relating to the same. However, 
on estimation of the parameters of the period of post partum 
amenorrhoea one may fit the theoretical distribution of 
the P.P.A again in the empirical data to get an idea of 
the different sorts of response error in the data of P.P.A. 

There are many other examples illustrating the 
indispensabiiities of stochastic modelling in the estimation 
of certain Demographic Paramters. For example, in a 
family planning action research programme it may be useful 
to estimate the proportion of female population at a time 
who are exposed to the risk of conception (fecundable 
period) and also what is the extent of the same population 
under infecundable condition (Biswas {1975)j for distribution 
of family limitation packages and to work out the cost 
of the same; which is impossible without appropriate stochastic 
modelling. Similar is the case of estimating the number 
of births averted under a given family planning action 
research programme andthe cost benefit analysis of the 
same. 


Finally, apart from exhibiting the applicability of 
stochastic modelling and almost indispensibility of the same 
in relation to certain problems it is worthwhile to mention 
the vital role of stochastic modelling in analysing incomplete 
demographic data based on foUow up studies. Either the 

reference period of follow up study is fixed in which case 
the proportion of observations having the characteristic under 
study is a random variable; or when follow up of the survey 
is terminated only on collecting a fixed desirable number 
of cases having a given characteristic the period of termination 
is itself random variable. Thus we have censored samples 

of type I and type II respectively. The analysis of incomplete 
data can only be done by appropriate stochastic modelling 
of censored distributions. 

In what follows the motivations of the models were 
thus intended to develop Mathematically oriented Demographic 
research and to support the Biomedical research concerning 
the biological factors affecting fertility, maternal and 

child health. While in India emergence of National problem 

of family limitation (because of rapid decline in the mortality 
rate with almost unchanging status of fertility condition), 
no doubt, gave great impetus In the development of probability 
models concerning human reproductive process; to evolve 

optimal strategies for the reduction of level of fertility 

on a natural scale but reallstrc Jiio^ befitting the assumption 


which the population really conforms remain a tough exercise 
because of basically two factors viz. (1) the lack of adequate 
data base to provide the necessary ingredients for appropriate 
development of the model, (2) Mathematical exercise is discouraged 
v/ith the replacement of probability models by simulation 
models presumably due to the ability of the latter type of 
modelling to take into account of several plausible factors 
which as an exercise with Probability models may still be 
a challenging problem. Nevertheless simulation models does 
not help much in defining a theoretical basis between the 
explanatory and' depending variables in the process. 

1 . 2 A historical survey of the stochastic fertility models 

and the motivation of this study . 

The first systematic study of modelling fertility analysis 
perhaps started from the paper of Gini (1924) which was presented 
at the International Mathematics Congress, Toranto where the 
notation of Fecundable was first introduced. In recent decades 
there has been considerable amount of interest generated in 
the analysis of pregnancies of cohorts of married women. 
Special emphasis of the pattern of live births, is one of the 
factors studied by Henry (1953). His work (1953,1957,1961) 
on models on Human Reproductive Process comprised of both 
discrete and continuous model postulating both homogeneity 
and heterogeneity of the population. 

A historical review of the models that they differ only 



in treatment of time as "discrete 


or "continuous 


many 


of the models the probability of conception is held constant 


and further assumption is made about all conceptions leading 


to live births which are associated with a fixed non-susceptible 


period covering the gestation period and the period of post-partum 


amenorrhoea. However, new type of more flexible and realistic 


probability models have also been evolved by introducing 


the notion of time when a conception is recorded. Emphasis 


is laid on the results for the probability of a recording at 


a specified time "t" in the same way as it is held in physical 


Muller Counter Model of 


The Human Reproductive Process is thus compared with 


Geiger Muller (G.M) Counter where every event occuring with 


an intensity is followed by a dead time (or blocked time) 


during which no further success is possible. The term "counter 


is paralysed" is used in this respect which is comparable 


to the non-susceptible period following a pregnancy in the 


The G.M Counter is classified 


human reproductive process 


according as the dead time is fixed or random variable. Type 


II G.M Counters are contrasted from type I in the sense that 


not only there cannot be any success during the dead time 


in type II (as in type I) but a result of arrival of a particle 
in G.M. Counter type II during the dead time further prolongs 
the dead time without registering the event during the blocked 


period 





For the analysis of human reproductive process, counter 
model I with fixed or random variable dead time appears to 
be more pertinent. It will be shown in the subsequent part 

of the theory that process (e.g. Dandekar’s model.etc.) are 
directly derivable from the existing results of G.M. Counter 
Models. Further, the renewal theory technique is amenable 



for the development of probability models on fertility and the 


human reproductive process. All these models can be classified 


The first deals with the waiting time 


into two 


distribution between two consecutive events (such as births 


while the second deals with the 


abortions 


probability distribution of the number of events during a fixed 


time 


In the first case, time is a random variable and the 


system conforms to ordinary or delayed Renewal Process, which 


In the second 


Non-Markovian in character 


case, the number of events during a fixed time is a random 


variable and conforms to a renewal counting process 


obvious to find out the former is continuous whereas the latter 
is discrete and one can directly be obtained fromjthe other. 

In India, the formulation of the stochastic model representing 


the human reproductive process was first done by Dandekar 
(1955). Dandekar assumed a one to one correspondence between 
pregnancy and live birth and a constant time independent Poisson 
intensity for births/pregnancies of all others. 
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A fixed non-susceptible period (dead time) was assumed 
following every pregnancy/birth in the pattern of Counter 
Model Type I. The probability distribution of the number 

of events in fixed marital duration (0,t3 has been obtained 
from the Modified Poisson Distribution as evolved by Dendekar. 
However, it may be seen that the distribution considered by 
Dendekar is nothing but a Poisson Renewal Process where 
every renewal is followed by a dead time of fixed measure. 

Thus Dandekar's distribution is otherwise obtainable as an 
example of Counter Model Type I with fixed dead time as 
considered by Takacs (1960), Dharmadhikari (1964), Biswas 
(1973). Further the difference equation techniques adopted 
by Basu (1955) provides alternative derivation of asymptotic 
pregnancy rate which is otherwise obtainable from Dandekar’s 

Model. 

Brass (1958) obtained Negative Binomial model from 
Dandekar's distribution ignoring the non-susceptible exposures 
following every birth and obtained a theoretical model for 

the probability distribution of births. However, the model 

did not give good fit in empirical data as infecundable exposure 
following every birth was ignored. The probability of zero 

number of births to mothers with completed reproductive span 
as obtainable from William Brass Model did not show good 
correspondence with that of the empirical frequency for the 
same' class because ^ of ' -(''Jp-ilesence ' oi con^i-:rablO' proportion;;,:- 
' of Biologically sterile' .'jmotheii'Jfbi'^ '-.ooiiplis ) ■ - , as - a ’ ! ^ , ' ■ 

' ■ 7 -;/ ^ 
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good proportion of mothers (or couples) who wanted to remain 
effectively childless without being biologically sterile. As 
William Brass Model made no provision for the same. Brass 
(1958) while attempting to improve the fit used a "zero truncated" 
Negative Binomial distribution which implied considerable loss 
of information. Hence the Brass Model needed further generalisation . 

Biswas (1973) made an attempt to generalise Brass 
Model by removing some of the point of deficiencies of William 
Brass Model is (i) ignoring the non-susceptible exposure following 
pregnancy, (ii) ignoring the biologically sterile couples which 
got mixed up with fecund couples in the population while 
accounting the distribution of live births. 

Dharmadhikari (1964) while attempting to generalise 
Dandekar's model, first cf all developed a model rased on 
Mixture of Exponential distribution with the origin shiftedjby 
K (the length of the non-susceptible exposure). Besides this 
generalisation of Dandekar's model modified into a mixture 
of two distributions, Dharmadhikari (1984) developed a generalised 
related time dependent process (Doub^ Stochastic Poisson Process) 
and considered the application of the same in further models. 

Singh (1963) independently obtained two different 
interrupted probability models on similar assumption with 
the problem of estimating the number of schools of fishes in 
a fishing ground. For example, the dead time "h" following 
each conception was comparable to the time for which further 
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scouting of fishes was temporarily abandoned following the 
sight of a school of fish and exactly "h" hours were spent 
for fishing the school. In both the models Singh assumed 
the presence of (1-^) of . the couples to be effectively sterile. 
While in the first model Singh assumed the monthly probability 
of conception to be constant, in the second model he assumed 



the same admitting a probability distribution 

following Potter and Parker. S.N. Singh (1964) 

mod.e-Lo ' 

thfi assiimotion of the piesence of a fixed 


; '■) further obtained 

coia/i is ' fe. 

proportion of effectively 

sterile couple as v.-ell as fixed non-susceptible exposure following 

every conception and obtained BAN estimators of the parameters 

of the model by using Central Limit Theorem. Modified Minimum 

method has been applied to estimate BAN estimators. 

As noted by Sheps and Menkin (1969). Singh’^ model is a modified 

restatement of Dandekar's model using a Neyman's Method 


(1949) 


Another aspect of S.N. Singh's work in this direction 
■s of analysing the building up appropriate probability 
for explaining the random changes causing variability 
nf rnnnentlons during a fixed fecundable period. 


A partial estimation of some 



more parameters of the model, 
of the parameters was made while assuming the knowledge 
relating to the rest. 

Perrin and Sheps (1964) made an explicit use of Renewal 
Theory for building of probability models on human reproductive 


process that allow a variety of pregnancy outcomes (live 


birth', still birth, foetal wastage, sponteneous or induced 



abortions) and general distribution for durations of pregnancy 
and the post-partum non-susceptible status. 

The states of the reproductive crtcess in Sheps and 


Perrin consist of 


Non-pregnant fecundable state 


Post-partum infecundable state associated with 
abortion or early foetal loss 


Post-partum Infecundable state associated with 
still birth. 


Post-partum infecundable state associated 
with live birth. 


The length of stay in each state is viewed as a random 


the waiting between two renewals of the same 


variable 


state and the number of renewals of each state during a 


fixed term (0,t] conforms to a Semi Markov pattern. To 
study the fertility process, it is, important to know the 
number of passages into : ft specified state, the fertility 
rate ’(or the pregnancy of time); /'and ^ ^ the 


pPi; 
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probability distribution of the renewal of 


different states 


of the process in a given period of time. 

The results obtained from these models have implication 
both for the analysis of data on the distribntlon of births 
to a group of women In a calendar period and for the 

planning and evaluation of family planning programmes. 
In this direction a family building model has been set up 



A very important 


and Radick (1969) 


Menkin 


model for studying the short term eltects or cnaug. 
reproductive behaviour t Sheps and Menkin (197X) ) on 

intrinsic birth rate and a model for conceptlve delays in 

a heterogeneous group have also come up. 

Basically Stochastic Models relating to fertility 

analysis can be broadly classified into two catesones; 

Viz. (1) Models Which predict the number of events (a discrete 
random variable like births. pregnancies. etc.) during 


a fixed period of time say 

(ii) Models which provide the waiting time distribution 

between two consecutive events or between any two events 

say ith and jth event (3 > i)- 

one interesting exercise relating to the models 

of the first type Is to consider the total number of conceptions 
or births during the entire marital span of a particular 
marriage cohort and compare the same with another proximate 
marriage cohort (differ lej by few years) having the same 
marital age while other sooio-economlo factors remaining 




more or less the same (Chapter-3). The motivation in such 
an exercise is to find out whether there is a change in the 
total fertility level of two marriage cohorts. Since the problem 
relates to basically two marriage cohorts where dates of 
marriage of the first and the second cohort are assumed 
to have taken place at two points of time with certain time 
lag; naturally the technique of Stochastic Point Progress 
was required to be adopted. The same type of problem has 
been solved by an entirely new technique known as "Palm 
Probability" which was hitherto used in "Queueing Processes" 

( Khintchine (1966) ) while facing the second type of problem 

on Stochastic Models on Fertility Analysis "Palm Probability" 
also happened to be useful to provide innovative solutions 
to some problems of dependent processes arising on account 
of the individual variability of the fecundability parameter 
following certain probability laws. Stochastically the variation 
of the fecundability parameter thus leads to a complete destruction 

of the renewal process, where each renewal corresponds to 

a conception according to the order of the renewal. 

The variability of the parameter leads to the result 
of interarrival time distribution between different order of 
conceptions which are correlated. Thus the waiting time 
distribution between 1st and 2nd order conceptions cannot 
be treated independently of the waiting time distribution 
of the 0th (date of effective marriage) and 1st order of 


conceptions. Hence it has been envisaged in our exercise 

to employ the technique of Palm Probability to obtain the 
conditional waiting time distribution between the first and 
the second (or higher) order of conceptions given that the 
first has taken place at a particular time say T=t. The 
theoretical development concerning this kind of exercise has 
been given in Ghapter-2 and Chapter-4 of this thesis. In 
the next place, a Multistate Markov Chain model has been 
developed to study effectively the efficacy of a particular 
sterilization policy motivated to reduce the birth rate in 
a population subject to the cultural acceptability of the 
same in a population (Chapter- j). A further extension 

of the problem of sterilization of mothers based on number 
of arising a family planning programme using the approach 
of Martingales is taken up in Chapter-6. 

In Chapter-7 of the thesis a partial solution 

t 

two other problems of related interest in the fertility analysis 
has also been encountered. The first comprises of an examination 
into the impact of certain family planning action research 
programme. Suppose that with the onset of a certain order 

of conception the level of fecundity of a population is reduced; 
then the problem comprises of observing the impact of the 
same reduction in the successive spacing between consecutive 
order of conceptions as well as to measure the reduction 
in the average number of conceptions following the implementation 


of the scheme 



The second problem _ comprises of re-examining the 
efficacy of the traditional post-partum abstin&nce as a social 
custom vjhile reducing the conception rate as an indirect 
measure. In fact, the period of abstinence is carried on 
mostly during the period when there is otherwise a post 
partum infecundable period due to lactational amenorrhoea. 
Therefore, the motivation in this study is to obtain the probability 
that the period of post - partum abstinence exceeding that 
of P.P.A. by a specific non-negative quantity. In fact it 
is worth-while to consider the difference of the distribution 
of the period of abstinence and that of post-partum amenorrhoea 
which may highlight the role of effectiveness of post-partum 
abstinence. 

The chapter-S of the thesis again employs the 

method of Palm Probability for obtaining a generalised probability 
model for measuring inter conception delays. 

The chapter-^ deals with a new application of 

survival analysis technique for obtaining inter live birth 
intervals by employing a Hazard model which can take into 
account of any type of fertility behaviour which one may 
theoretically assume. 

The study concludes with a new investigation which 
may highlight the interrelation between two useful variables 
viz. the period of post-partum amenorrhoea and the lactation 
period; which are useful items of consideration for averting 
birth without use of mechanical methods. This has been 
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done by a new application of renewal models which 
has been, hitherto, used in survival analysis, using 
Freund's Bivariate exponential model attempt has been 
made to construct the regression of Post Partum amenorrhoea 
period (P.P.A) on lactation period (when P.P.A. 

> lactation period) for prediction of P.P.A; as well 
as regression of lactation on P.P.A. (when lactation 
period > P.P.A). 



The practical importance of 

Freund's 

model 

in 

this 

application is because of the 

assumption 

that 

with 

the 

ceasation of lactation the hazard rate 

for 

the 

discontinuance of lactation is also 

increased 

given 

that 


the P.P.A has expired. 

1.3 A Survey of the Outline of some of the Methodologies 
employed . 

1. Palm Probability is defined as the conditional 
probability of a specified number of events in a time 
interval given that an event has happened at the beginning 
of the interval. Cox and Isham (1980) describe Palm 
Probability as follows 
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For u < V, let N(u,v) be a random variable giving 

the number of events occuring in (u,v) and X^be tte sequence 

th 

of the intervals between i successive events (i=l,2,3,.....) 



in a process starting from an arbitrary point 


Consider the survival function 


P(X > x), where X is a r.v. representing 

the waiting time for first renewal 


to orderliness which irr.olies 


the infinitesimal interval) which is the limiting 


renewal in 


that, given an event occurs immediately before 


probability 


the next event of 


occurs 


instant (i.e. the system will survive more than a period 


Now, by stationarity condition 

? [N Co,x) - 00 c- S.o') >ol 

« f CN Co, '=9 o] - p c Ki (- s; X) ~-ol 

>= ? L N Oc) - 0 ] - ? ) -o] 

p [Ni Co,'x) = o I N ^>0) > o] 
p r M TOiX') =.o A Nl C*" ^ ^3 




P =0 1 M C-S,o>ol 8'' PtM CS)>o] 

g-' j-p (x-) =c;i-T[NCx + 5') = c13 

g-' CPCNCx.+ S)=o]-PC^tC-)=°^]• 

£rm. 5'' p [N CS) >o] 0-5) 


as the 


occurrence parameter ^ of the process which we assume 


as finite. 

Fa9tiKex.? dt£ndt*.rv^ 


-V ( 1 . n - r, -x'', K -- - ’ 

, 1\ 'X'V -• J " ' \ 


i u 


' 4'^ 


as 


the distribution of N (x). then in the limit as 


^ ii.u 'p [ N Co,:c.) =ol N C-S,0')> 0'] , v ) 

04* ' ' 

No"' 'P[N(o.>5=ol N C-S,o')>63 ^ ■P (S>>ol 
== K Cx)r' ? CnCS)> 0] -fn<n«-0-+'> . 

X 
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Since the process is orderly, if k ^ 0, as S — > 0+ 

= ?EN(-f. O)=o, MCo,=<-:) = k] 

4- PUN C-S,o') = \, M Co,x> = k-i’] + oCO 
= U C^) — P CM (- o) > q, M (o^5c.) =Kl 

+ ? CN C- S', o) > O; N (o,x> = krO-^oC^) 


SO that 


f> (x-t-g') - b c-^) = - p [N C- £jO»c , H fo,x.)-lsl 

^ p'^rx g,o-) > c/n Cc, X) = k-il l-OfSj 

- r K,' ■' /I r- !”iM ^ < M 

+ P [N C-S,o)>ol. PtN (:o,x>= k-i| N(-S>0'»6 ]->-o(;s3 

+ s-'pLH(-f,o)>o].?i:N(o,x>k-i I N 

__£■'[.? [NC-SjO')>o3.T'Ct'i(o,x>klN(-£o)>o] 

- P[N(-S,o3>6].?CN(:o,*) = k-iJ n(-S,o>cO 1 - 0 ( 1 ) 



denotes the derivative 


where D, 


The equation il.5) links the distribution of the interval 
between successive events with survivor function 


of forward recurrence time with survivor function 


pCx) (where P is the probability that starting from 

0 0 ^ 
an arbitrary time instant, there are no events in the following 


interval of length x 


as a limiting probabilit 


Palm Probabilit 


e general results connect 


Integral Equation ) 


distributions of events conditional on a point at the origin 
with those where the origin is an arbitrary instant may ue 
obtained. Assuming that the process is completely stationary 

and has a finite occurrence parameter and is orderly 

so that i) is equal to the rate of the process. 

Then tor each x > 0. the palm distribution is a discrete 


distribution defined by 


In a careful mathematical development the existence 
i) and ‘ir lx), and more generally of other limiting 
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probabilities of events B given NC-d.O) ^0 of the form 

itTV p [B 1 N (-£, ®) > o] (i.v) 

has to be proved. 

Let the probability measure H be defined for events 

B on those processes which have a point at the origin. The 
measure H can then be shown to satisfy 

II(3)= liw Tt'E’lN 

^ O-h 

for a wide class of events B. The measure H is called 

the palm measure of the process. In equation (l.S) the distribution 

of the interval measured from an arbitrary time 

instant to the next point of the processes linked to that of 

the interval bet\veen successive points. In the same way 

the functions ^ (x) defining the palm distributions given 

■ k 

in (1.6), which specify the distribution of the number of 

events in an interval of length x given a point at the 

origin, can be connected with the functions P (x) which 

■ k 

give the distribution of the number of events in an interval 
of length X starts at an arbitrary origin. These connecting 
equations are known as Paltn-Khin£chlne equations and may 


be derived as follows : 



Hence 


lienee- oA q^ 

dk C |\Cx)] = §-f p r Ni c- r, o:) > cl • C “ 

= - [pj. DO ~ '^-1 f^Oyn (/.3) 0.1 

where D denotes the right hand derivative. Th( 
equation for k = 0 has already been derived and is 


The integral forms 


It follows from (1.9) and (1.10) that 


::24;; 


Therefore, the probability of having exactly k events 
in 10, x) starting from an event at 0, can be obtained by 
differentiating the probability of getting not more than 
K events in (0,x) where 0 is an arbitrary time instant, 
Alternatively^Cl.ll) arid (1.12), 

0^) 0c)+-- - ■ ' - ■■ + (pO ^ 

=z 1- du- 

so that the probability of getting not more than K events 
in (0,xl, where 0 is an arbitrary instant can be obtained 
by integrating the probability of exactly' ,y events in the 

interval when there is an event at 0. In addition the 

right hand side of (1.13) is equal to 

no 


0 - ' 3 ) 


J 7r^(c4)<ia- = Aj TU^Cu-') cite 




CP 


j:? ^oO) + |>(oO])f - - ■ ■ b I — “Pv J TI^ {^(^3 


oo 


J -TV (ft) d 


u =■ i 




fkAce. + = ? Cm Cxi 4 J 

(1.14) may be justified by the argument that if 0 i. 

an arbitrary time instant and there are no more than 

u PVPnts in (0,x) then there must exist an event with 


y p? 0, such that there are exactly 
Since the process is orderly, the 


coordinate \ 


of an event in 


The equations (1.9). (1.10). (l.lD and (1.13) can 

be summarised by using probability generating functions. 


For if it is defined 


^ ■■ ■■ 

arbitrary origin while G refers 




to the situation given a point at origin 


hand side denotes the cumulative distribution 


where the 


obtained by k-fold convolution 


That is 


Birth and Death Process is a process which permits a population 


as well as to decline. These are more relevant to 


biological population in which both births and deaths occur 


The birth and death process developed here is applied 


to some problems of policy making in sterilization in chapters 


Let X (t) - Size of the population at time t 
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with the initial size XtO) = k 

.0 

and P (t) =0 where X (t) = k. 

k ' 

Given Xtt) = k we assume that the probability of exactly 
one birth in the interval ( 4; o(^!) 

The probability of exactly one death is 

o(£) 

The projability of more than one change is o(S)- 

Therefore, the probability of no change in (t.t +?t-_) 
is 

Consequently, the probability , a. so.e 

may be given by the Kolmogorov equation 

•)-oCO 

This given us the corresponding system of differential equation 



(i) 

(ii) 



:; 28 :: 


©1= Ct) 

It' ' !-■ t ^ - 'K, ■' '' 


4 k ft) = - i)\^ft) +-A ^ 


Under initial condition p i 

bQo) = i ^^Co>='=> k=kc- 

( 1 . 19 ° completely determines the probability distribution 

Appropriate assumption may be made regarding the 

function it) and It) to obtain stochastic process 

k . 

corresponding to empirical phenomena. 

A particular case : linear growth 




O-O-o) 




-{(kri) ■ A ■ ft) • 


(l-O-l) 


Employing the method of probability generating function 
ip.g.f) to solve the same j 







From (1.21) it follows 


Viith the initial condition t 


The auxiliary equations are 


rtial fractions to rewrite 


the first auxilliary equation as 




The second auxilliary equation 


Q = Ccmcha-rk 



solution 


particular 


given by 


Therefore the general solution of 11.23) is 


where (h, is an arbitrary differentiable function. 

Using the initial condition, it can be seen at t 0 

k 






at least for all s >vith /sj 4:^ 1. 

Hence for all 9 such that 


(h [ 0 ^ = ^ 

1 S 

■ i 

\ + 0 A 


(l,30) 


Putting 




K 


The p.g.f. can be rewritten as 

r o({0 H ’ 


\ ; 0 




1 


,'3f0- 


(/ 3 s) 


Expanding the p.g.f- G (s:t), 

.X , 




-f. f Xlf) -- 

L_i 








f’K ^ Z 


4 


c((t) 
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The use of Martmsales based on a density dependent birth 
and death process is done in chapter 6 on sterilisation policy. 


A Stochastic Process X 


n = is 


n 


a martingale if for all n = 0. 1.2.3 

. ,0 e{ ' 


oo 


hold 


(-0 ^ i I 

A more general definition *. 

Let X : n = 0,1, 2, 3 

n 


■]- 


^ (rU) 



where 


exists for some 


distributed random 


identically 


are independent 


then 


is a Martingale 


Q„h Martlngal- Martingale 


is G&llod ^ Supsr 


be stochastic processes 


and 


( Y i n - 0 , 1 , 2 , , • 

: n 

We say \ 


) are stochastic processes, 
is a Martingale ;vith respect to 


if for 

all n = 0,1. 2, 3 


( i) 

e ( 1 K^O < 


(iO 


^ Y n n 

One 


Y ) as the 

may imagine i x » ^ • . . • • 
0 1 

n 




history upto the stage n. 
fjald's Martingale : 
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Martingale with respect to (Y ) if for all n 


(0 

£(x;) > - « 

X - i'rf [ X/°i 

in) 

Y, 


(iu) 


7ipy) of Yfl ^ Y » Vt 


Similarly (X^ ) is called a Sub Martingale with respect to 


(i) £C Xn j <■ ^ - Suj) (x, <9} 

(ii) B l ^ 

(ii‘' {5 a ■iuv.cfion YoaV, .? 


Definition of Markov time (or St opping time) ; 

A random variable T associated with a Martingale 
(X ^ is called a Markov time with respect to (Y if T 

takes the value 0.1. 2..-. and if for all n = 0.1.2,.... 
the event T = n is determined by (Yq .Y ^ .....Y^). 

By the term ’determined' implied that the indicator 

function of the event (T = n) can be written as a function 

of Yjj . Yj .Y2 
Optional Sampling Theorem : 

Suppose (X is a Martingale and T is a Markov 

time with respect to (Y^) then 


:;34: ; 



rK-^h { i' 


interctiange 


can justify 


whenever 


Thus 


following 


limit 


holds 


nptinnal Stopoing_lllgo£ga 


be a Markov time 


^ .eprese..s the Indicator iunotion tor n. 

where X h V 

. ThPorv is employed to evaluate 
in Chapter 7 Count^WlSfill ta P . „t 

^ . 40= for the reduction of 

. inrtlrect strategies tor 
the direct and mdirect 


conception 


Mnapi I is a registering 
Counter Model t 

of a radio active 

detecting the presence of 


Muller 


mechanism 
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because of inertia, 
the impulses- More 
arrives at a fixed 
The registration 
and impulses during 


counter 


some 


counter 


an impulse 


sters the impulse 
of length 'Hi 
he registered by 


time L ^ 
causes < 

the 

general, 
of the 
counter 

Hm-jever in the simpler cases 

"" . a de'aa ti.e do not cause any dead 

which arrive during a - bv a registered Impulse. 

, at each dead time is caused by a registe 

,bis is caued a Oeiger-MuUer Counter ct type i. 

in a counter ct type U (Smith i«) -h am 

■ ise causes a dead time sc that arrivals during a pen 

dead time. A counter model 

description of certain physical processes. 

Biological, social and Industrial 


the counter 


time 


termination 


> after 
will be 


first impulse to 


registered by the 
of length say 
the impulses 


time 


time 


andthis again 


time so 


time prolong further the 


although basically u 
can be applied many 


fixed dead time x 


mndp.l of type 


Counter 


I with fixed 


„e Illustrate the counter model type 

dead time using the foUovilng problem : 

* conception takes place with intensity 

re the condition that ever, conception Is followed h 
exposure w (fbced dead lime) during whi. 

. te.no Then to obtain the probabi 

conception takes place- Then to ou 

’'of' the' nitmlJer: 


Let 







(Poisson 


(waiting time of conceptions) with conception rate 
intensity) having negative exponential density function 


-tLi^ 


Probability of not more than n 
(conceptions) up to time t 


of the random variable 


where 


the corresponding survival 


function 


Denoting L(.) as the Laplace transJorn, 
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( 


L(F 




(tj = _L Llf 

5- 


(t)) where 
»th 


n+i 


n+1 ~ n-l 

(.) is the interval density function of (n+l3‘'“ renewal 

(order of conceptions). 

I. (t> f 


nfi 


where 


is the density function of 


y f 

'i } ^ 


is the density 


function of 

(\ - r^. 

J, 


Ci n<j 

* 

Stfxnds 

4or 

. f 

c<5> *1 vol , 



111 

5 

•fold 

CO nvoiutio y) ® 


L 

<*') L 



1(1^ (t]) 


A 


VI J L ( ^ (t) ) - J 


A t S 

oa 


-St 


TT 


X 


-ITS 


A ^ 


S4 A 


C df 


(hH6) 




A 


nf I 




(5-t A) 


^h^ I 






Putting (1.46) and (1.47) In (1.45) 




Laplace Transform 
■f - ri 


Also using the result 


Putting 


Further on substitution 



Comparing (1.51) with (1.49) 


where 


refers to the greatest integer contained 


Counter Model Type I with Random Variable Dead Time 


We illustrate the concepr of counter model type I with 


the following problem 


A particle arrives at t=0 and locks 


Problem 


With the registration 


the counter for a dead time of duration Y. 


of the particle the counter is blocked for a time of length 


The next particle to be registered is that of the 


first arrival once the counter is fixed. The process is repeated 


where the successive locking times denoted by Y 


are assumed to be independent with a common distribution 


and independent of the arrival process. We obtain the waiting 


inter-arrival of the 


time distribution of the 




Where 


is the dead time following the 


is the residual life time before the counter is locked 


Where Y is the dead time and Y be the residual 


life time before the counter is locked 


Denoting 


the density functions 


of the dead time and the residual free time (before the system 


is again being locked) have 


Counter Model Type II with Fixed Dead Time 


Here the locking mechanism is more complicated 


As before, an incoming signal is registered if and only if 


it arrives when the counter is free 


counter 


only recorded particles Induced the counter to lock 


For Type II Counter every arriving signal can prolong 
the dead time in the counter, the associated locking times 



: : 41 ; : 


being added concurrently 


a 


hfrlvJ iy^tystr<^'t\oy\) 


Cl ^ XA.f' t. 


T / T <^A are the registration times 

Li0X 1 I • 2. 

^ onH thp nth renewal and arrivals take 

of the first, second and the ntn reaewci 

\ Thp Isnstli of ttiB dBSci tirns 

place with Poisson rate A 

following every registration is iT- 

c /V /V t are the independent identically 
We assume^ 1 ^ - I j 

distributed random variables. 

T -AX 

p[^ ^ e. 

Pl'i'v,- V, ^ 


where the structure of F(x) is to be obtaioed. 

<p(^) = I = llk-t)] 

Let uit) = dUU! be the renewal density of the process and 

dt 

U(t) be the renewal function. 



the Renewal Density is given 


Then the 


oceed as follows 


To obtain 




that initially the counter is free 
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df 

0 . d 'tr- ' 

- Am At-^ U-^^JmAr 

o di: 


-Ar 


L^U(t-i) 


ir oo 

; ^ 

, - x^A+Sjr- _(;t+s) 7 r 

> e __ ■ A €, _4 _> 

~ ^ " * A-vS A + ^ 


- - 


u (t) = dU M 

dt 


'■ -*1. ■ 

L I 


X + 5 


A + s 
'Ar 


L ( A t 


) -ye 

f^ir ,J__ J_]V 

ktS/ J 


^ S At 5 


A e‘ 






ATT,. .-Xf*--"") 

iy>T^ 


A e -iri--^ 


■ 



by the second shifting property of Laplace Transform 


(MS)jr 


Counter Model Type II with Random Variable Dead Time 


This counter process is quite difficult to analyse in 
general form. However, the results are known only for process 
with Poisson inputs having arrival rate ^ . 

Let pit) be the probability that the counter is free 
at a time t/i and registration is possible. 

It is to show that 


I 




Gty) represents the cumulative density function of 


where 


the dead time distribution 


Given 


occurences 


the distribution of occurence times is the 


interval 


independent random variables taken 


from a uniform distribution in 


The counter is free at time T 


Proof 


if all dead periods engendered by these signals have been 


terminanted before t 


P (dead time commencing in y will 
end before the time t). 


Since the locking times are assumed to be independent 


and also independent of the arrival process 


We have 




But the number of signals arriving during the interval 


has a Poisson distribution with mean 


From the Idw of total probability 


^ p(t) is the probability 


Counting with our assumption 


density 






dU(t) 


the renev^al function in 


where 


phS t 


U(ti- pfr )Jt 


Putting (1.71) in (1-72) 


; _AlO-&CYjMx 


the result 


iVi'i 




CHAPTEB - 2 


OH THE APPLICATION OF PALM PROBABILITY FOR OBTAINING 
INTER-ARRIVAL TIME DISTRIBDTION IN WEIGHTED POISSON PROCESS 

2.1 The distribution of = I • n = */ 

where T^ represents the time for the n^^ renewal in a 
Poisson process with parameter A is given by the density 

function. 

%{ trl X) .* A 6 i •» <- i:< j A > o 

However, if A , instead of being a constant assumes a 

probability distribution with density 

^ k — a A U-\ 

then the unconditional distribution of time till the first renewal 
lor arrival) is given by 

( f(tl»y(A)dA — (2-3) 

, th 

But the waiting time distribution for the n arrival given 
that the first took place at T*t is not certainly the n-fold 

convolution of (2.3). Because of weighting, the renewal 
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structure of the Poisson process is completely destroyed leading 
to interarrival distributions having infinitely divisible structure 
but with dependent increments. Denoting by X(t) the number 

of events in (0,tj from a Poisson process weighted by a Gamma 
distribution given by (2.2) it follows that y t ^ s 

~ f [>((s) Kit) - y(^)Ix] -[eiMjsti-s) 

3 s(t~s) variX) >o {%-H] 

This makes the problem of obtaining the probability 

distribution of the interarrival time in a compound Poisson 
process (or any such compound dependent process) somewhat 
compiicated. On the other hand such interarrival time 

distributions are often considered as very useful for practical 
purposes; say while obtaining the distribution of interpregnancy 
(or interbirth) intervals or the waiting time distribution 
between two morbidity spells or the interarrival time distribution 
between two accidents. In all the cases the hazard rate 
j^(t) (or the intensity) even for the same t varies between 


individuals. 
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Cox and Isham 11980) have dealt v;ith a series of such 
problems on dependent process classifying the problems according 
to the nature of dependence. This chapter is devoted to exhibiting 

an application of Palm probability to obtain the interarrival 

' til 

distribution between the first and the r arrival (r = 2,3,.-....} 
for a dependent process. 

2.2 The Result ; 


Defining Palm probability (t) 


probability of 


events (r 


as the conditional 
0,1, 2. 3 .) in (0,t;] 


given that an event has occurred at time T=0, Khinohin {I960) 


giwes 




oin A 


where 


v; (t) 


(t) 


if] df 


(X-() 


represents the unconditional probability 


of r events (arrivals) in (O.t] following a weighted 

Poisson proc<iss ( weighted by Gamma distribution as in (2.2) 






i., 1.- 


1 - V. Ct) 






Is the intensity of the process (Khintohine, (1960^. 









Then 




Using 12.6) and (2.7) we have 


which on differentiation gives 


an event has occurred at 


is the random time of the first event 


Where 


following the occurrence of an event od. T=0 


Hence 


is the conditional density of the waiting 


where 


time of the second event given that the first has occurred 


different from the unconditonal 


waiting lime density of the first arrival given in (2,3) 


r=l and differentiating 


Procoae(ing in this way, putting 



both sides of (2.5) we have 


Then 


)is the conditional cumulative distribution 


where F 


function Ic.d.f) of the waiting lime distribution of the second 


given 


Then 


K(K-t a 




which is the interarrival density function between the first 


and the third arrival given that the first arrival took place 
at Ts=0 in a compound Poisson process. Proceeding precisely 


in the same way we have 



t(trth) 


when, we get 


Using 


relation 


(JL f f 
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I I I 

t (K + l ■'•(Kti^) 


(►I-')! fa + t) 


*•* (2.H} 


which provides the distribution of the time between the first 
til 

and the (n+1)^ ■ arrival given that the first arrival occurred 

at T=0 for a compound Poisson process weighted by a Gamma 
distribution. 

The treatment for other dependent processes is precisely 
the same as illustrated in the foregoing example. 

Remarks : 

th 

If IT ^ - T _ ^ ) is the waiting time for the i 

order of conception (i = 2,3,.....) or (T ^ ~ T. _ ^ ) may be 

called ‘Interconception' interval and the conception rate 
^ (even taken independent of i) varies from individual 

to individual conforming to some probability distribution say 
Gamma distribution, then because of waiting ^ 

(i) Renewal intervals “ ^i-i ^ ^ i = 2,3....... will 

cease to become i.i.d.r.v.® (or i.d.r.v.s.) contrary to the 

traditional assumption. 

(ii) The number of renewals even in two non overlapping 
intervals will be correlated. The extent of correlation because 
of the weighting of the process. The renewal structure Is 
completely destroyed; leading to the process to conform to 
infinitely divisible distribution with dependent increments. 


we have 



Correlation between i and 


Let us assume the probability distribution of 


We have 


E [ E Tj 1 A)]- e ( E (T| lA) l(E(gfl] lA^) 


CM t-TC ,TJ] 

ECt: Tj 1 a) E 


cfTllAj£fT:_jlA)l 


being two non overlapping intervals 


£Crr l^) 


(k-l)(K-i) 

£ ( r - T : ) - 

CovfTrTi)^ , 

c.v(r;-,lj UCova;,-rj tTj 



CHAPTER-3 


On the Comparison of Cohort Fertilities 


by Palm Probabilistic Technique 


3.0 In this chapter a suitable methodology is evolved 


to predict and then to compare the future fertility performance 
(vis-a-vis the total fertility status). Given the time of 


respectively 


occurrence of first birth as t 


of the two cohorts, whose age of marriage differs by a 
time lag say, < it is proposed to analyse the level of fertility 
changes between two cohorts by the average number of births 
to these cohorts of women given that the first birth to the 


and the first birth 


first cohort of women took place at t 


of the second cohort of women took place at ^ 

leaving the residual fertility span (T-t j ) & (T-t 2 ) respectively 

for the two classes of women who complete their fertility 


The purpose is to see whether postponement of the 
age at marriage is effective in reducing the fertility status 
by using the technique of Palm Probability. 



Palm ProbabiUt 


. ,vna involving pL'edicUon of the number 

Stochastic modelling involving h 

. fn tV or the prediction 

, events during a fixed period Of time fO.O ^ 

, rhe waiting time distribution between 
Qon anv two events say the 

events or betw 

event can be so conditional probability 

as -a'- „ .nat the event has occurred 

. oomber of events given that the 


of the inter 


consecutive intervals 
nr.:n in the follow in 


, rn s> be so small such that at 
let the interval (0,s) be ^ 

, of the Shoch with compound Poisson inputs Cl.c 

tor family of Gamma 

. t X weighted by a two parameter fam 
Poisson input teoundablU.y between 

K, and let US assume that there is 
" Js; without any loss of .nerelity. 

, 3-^0. the point Of errlval denoted by XUn PI.- 

rh llmiung posiuoh wUl be located In the beglnafn* 

. ,30 «hat IS th. probability of having 

nf the interval U.ti. 


Further 


distribution 

individuals) 


arrival m 
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a finite number of arrivals of shocks, say k number of 
arrivals (k=0, 1. 2, . . . . . ) in the interval (s,t) of any finite 
length given an arrival at the beginning of s? What is 
the waiting time distribution of next arrival or the nth arrival 
in sequence following that event (i.e. an event has occurred 


at s), 


These limiting conditional probabilities are known 


as Palm Probabilities. 

Let (t) be the conditional probability of k 

number of births in (0,t) given that a birth has occurred 
at t=0 (i^j^(t) is a palm probability measure) 

and V (t) be the unconditional probability of 


k births in iO.t); k = 0,1,2, 


Then Palm's integral equations connecting V (t) 
and (t) are as follows (Khintchine (1960) 

V,W= f K=o,U,..- 

n O ■ ■ 


. ■ V ■ 

v„(t) X l,(ti [ ^ (r )cir •■■13.^) 


where h (t) represents the intensity of the process that 
the hazard rate of having a birth at the age t. 

Since the hazard rate h (t) is a decreasing function J: 

/ ’ : . i ^ ' 









woman with fecundability level 


as the conditional hazard rate 


Assuming the variation of fecundability between women to 
women conforming to a Gamma distribution with two parameters 


have the unconditional hazard rate 


where t ts age of the 


Also 



Prob(no birth upto time 


r ii(ru)dr. 


where 


showing the stationarity 


which is independent 


Palm’s integral equation 


differentiating 


(a+A(t)l 


where A 






■ 



conforms to a differential difference equation 




Since 


given by n 


we have 




We have 


Again by Palm's integial equation 






U)kc/^ 






Precisely in a similar way 


('K+2) d 


xpected number of births in the residual fertility span 


that the first birth has occurred at t=t 


the first cohort is 


a, t A(T-t,)) 



and the expected number of births in the residual fertility 


span, given that the first birth has occurred at t=t 


the second cohort is 


where (a 


the Gamma 


ion corresponding to the cohorts I ana II respective! 


distrib 


Therefore for a change of age at marriage by 


there is change in the average level of births 


indicates fall in the overall fertility level 


indicates increase in the fertility level”, and 


indicates the stationarity in the fertility level 


Estimation of the Parameters of the Model 



cohort 


have 


birth at t 


second cohort 


) I there is a 


number of births for the first cohort 


Then the variance 


Var(XlT-t 


( t ■ \ 

Arr-t,) 


Acr-%) 






the second cohort 





the method of moments 


and T-t, will enable us to estimate a 


approximations 


Similarly by the same procedure } and 13.22) for 

given S and T~t - the estimates of a- and k 2 will be obtained 


by successive iterations 


On substitution of the estimates 


(T-ty) and tT-t„) one can examine 


the positiveness (or negativeness) of the expression in (3.20) 
to ascertain the increasing or decreasing trend of fertility 
as observed in two cohorts who differ by years at age at 
marriage. The inherent difference in the fertility, apart 


is 


from the availability of total marital exposure, if any, 
expressed in the estimated values of the parameters (a 
i = 1,2, for two cohorts. 

A Numerical Illustration 

Let 

E,(X(T-tj) I (t^) = f^(a^,k,,S.T-tj)=iiO-'3-64,4.48, 5.50 

Varj(X(T-tj)Jt^) = ^ =10.49,9.65, 10.69 

and = 45-16 = 29 years 

and 

E2(X(T-t2)lt2) = ^ = 3.09,4.40.4.62 

and 

Var2(X{T-t2)ft2) = 14 ( 32 , k 2 , ~ = 9.43.15.20,7.90 

and T-t 2 = 45-18 = 27 years. 

Solving (A) and (B) for the given set of hypothetical values 
as given above by starting with an approximate initial solution 
of a j = 3° and k ^ =kj. we get by successive iteration 

■ ■ ■ ■ 

c^^c 0 ,'i 5 ^ 0 ,^O / Kj r l.o y 1.5 j 2<^ 

Similarly starting with an approximate initial solution of 
a - a? and k, = k? we solve (C) and ID) by successive 

iterations and get 



Substituting the estimates 


ive get three values of 


0.48 respectively 


which shoves that in the 


hypothetical exercise the first cohort 


has the fertility status 


not significantly different than that 


of the second one 


A note on the nature of the Palm Probability Distribution 


of the number of births during fixed marital exposure 

II I i 


variance 


comparison 


mean 


Palm Probability distribution given in (3.173 shows that 


Variance = Mean + c. (Mean) - (mean) 


where 


(mean) 


relation so that 


For the non-negativity of the variance it Is necessary 
I for Mean 1. 



Variance = Mean (1 + (c . -(mean)) 


Variance = Mean (1 +(mean (0-1)) 


In particular it may be noted that if mean is greater than 


Variance •< Mean; or 0 


Variance 


Mean 


0 < V- Variance ^ Mean. For 6 = 2, Variance = Mean 

= (Mean)? 


Hence the model is flexible so that by suitable choice of the 


the mean/variance ratio can be made in the same 






GHAPTER-4 


ON THE APPLICATION OF PALM PROBABILITY FOR OBTAINING THE 


WAITING TIME DISTRIBUTION BETWEEN THE FIRST AND THE HIGHER 


ORDER OF CONCEPTIONS 


Singh (1964) has modelled a probability distribution 


on the time of first birth. Assuming that a conception taking 


place during the fecundable period following the marriage 


the waiting time distribution for the 


with hazard rate 


first conception is given by 


On further assumption, that the fecundability parameter 


varies from individual to individual (even with the same 


parity group) following a probability density given by 


the waiting time distribution for the first conception following 
the effective marriage is modelled by Singh (1964) as 


However, the problem of obtaining the waiting time distribution 
between the first and the wcond conception or between the 


9 



This makes the problem of obtaining the probability 


distribution of the inter-arrival in a compound Poisson process 


(or any such compouna dependent process) more complicated. 
On the other hand, this is the type of distribution we require 


while obtaining the waiting time distribution between the first 


and the second conceptions. Using the concept of Palm Probability 
(Khinchine (I960)). Biswas and Pachal (1983) have however 


obtained the probability distribution of the 1st and the 
nth conception (n=2, 3,.....) given that the first conception 


has occured at time T=0 
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first and the nth conception (n=2,3, ..) cannot be obtained 

as a direct generalisation of the result obtained by Singh 
(1964); because on account of the weighting of the Poisson 
Process by a Gamma distribution the renewal structure of 
the process is completely destroyed leading to inter-arrival 
distributions having infinitely divisible structures with dependent 
increments. (Biswas and Pachal (1983)). 

Denoting by X(t) the number of compound Poisson events 
(weighted by a Gamma distribution given by (4.2) it follows 
that for all t s 

Gov (X(s), X(t) - (s)) 

= (X(s) (X(t) - X(s) lA ) - CE U)J^S(t-S) 
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which is based on a similar assumption of having a time 
independent intensity of conception and assuming to conform 
to a distribution as in (3.2). The results have been proved 
in earlier chapter (vide Chapter two). 

For n=l i.e the probability of the waiting time for 
second conception subje t to the first conception taking piace 
at T=0, is therefore, given by 


fl(t 1 


) 


k+1 

a (k+1) 
k+2 
(a+t ) 


. ... (4.6) 


it may be noted that the probability distribution of 
the waiting tine between the elfeotive marriage and the 

first complete conception (given by (4.3) is obtainable from 
the corresponding distribution of the waiting time between 
the first and the second conception given by (4.6) replacing 
the parameter K+1 to K and vice versa. 

However there are several ways in which the model 
(4.5) can be generalised. As in Chapter 3 if we take 
h(t\A) = , &>0 instead of h(t 1 X ) = X based 

on Which (4.5) is developed, we have by the same notation 


as in Chapter 3. 


where 


G(z,t) - Probability generating function of ^ 
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the conditional waiting time 
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+ (k^-2)(k+3) fZMll ^2 


+ (k^2)(k+3) 


z(l+z+z +z + 


f i+(k+l)A(t)z + ( k+1 ) ( k+2 ) ! zA^)_ 

I "T ^ n '<■ a ^ ^ 


+ (k+1) (k+2) 


1+ (k+2)zA(t) 


+ z(k+l) 


+ (k+2) (k+3) 


+ (k+2)(k+3) 


(k+l}(k+2) fMlli 1 2 


+ fk+l)(k+2) 


l+(k+2) zA( t). 


z(k+l) 


(k+2 )(k+3) 


V. nf z"on both sides, we get 

Equating the coefficients oi z 
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is the survival function corresponding 
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Numerical Illustration 


^ as obtained 


in Chapter 3 viz 


the estimates of expected time from Marriage to first birth, 
first to second birth and second to third birth (assuming 
one to one correspondence between conception and birth) 
with the further adjustment for decreasing fertility level 
with age so the parity level of the mother are given 


as follows 



4.5 ^ 



Conclusion 


A further modification of the results when marital 
re T is considered finite is given by Srestha and 


Biswas (1985) 


Here the density function of the waiting time distribution 
from marriage to first conception is given by 


of conception under 


The monthly probability 


given by 


(a+j+k'T) 


In particular for j 



However, the result for advanced parity groups viz. (1-2) 


and (2-3) are still unavailable. These results if available 


could have really highlighted how much additional precision 


is again when the assumption of infinite marital exposure 


For the first order of conception the gain 


is waived 


in precision is considered to be negligible. As such the 


exercise has not much motivation 






CHAPTER- J 


A MULTISTATE MARKOV CH AIN _MODE^FOR_EyMiMMigg 
A STERILIZATI^O^POUCY 


sterilization Policy for the reduction ot toe 

rate in a population is often adopted on the basis of the 

number of surviving children. Generally, a sterilization 

* cftrar n rminle has a desirable number 
policy is implemented after a couple 

mina a certain hypothetical survival 

of children on assuming a cerrain yn 

rates operating in future. However, this method has certain 

degree of arbitrariness : viz. the long term effect of mortality 
in deciding the number of surviving chiidreu may not 

be properiv accounted for. ..n improved techuigue of assessing 

the long tern Impac. of such sterilization policies is therefore 
iu consider the Joint effect of f-ermitV 

The present exerctse has been made on the basis of the 
same rationale. We have developed a Multlstota Marhev 

Chain model corresponding to varying fertility and mortality 
intensities at different levels (States, of surviving children 
at the couple. The asymptotic probabilities of having a 
V. nf rhiiriren have been worked out. 

fixed number of chuaren uav 


novftinpment of the Model 


Notations 


n number of children 
(n»0, 1,2, 3, 


ProbablUty of 
unto a time t. 
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Probability of having (n+1) children 
in (t+ ) given that n nuniber 


of children are surviving 


time t- (n=0, 1, •<!, 0 .■•••••• 1 ■ 

Probability of having (n-1) children 
(t+ A ) given that n number 


) ^ h +o( A ) 


surviving 


children are 


time 


A possible upper bound of 


Number of surviving children 
at the time of sterilization based 
on certain policy. 1 $ 


the Kolmogorov equations 


Then we have 


) ( 1 - A ^ A -h 0(A) } 

[ 1 - 


(t -d- A) 
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and 


( 6.1 ) =4> p' C-t ) r - A. (-6 1 t fj ( t ) 


Now 


/ 

(5.1) 

/ 

'( ?• 2) 


P„^., 


I 




and 


(5-3)=$. Pm-k+i^’* .! 


staadv State Solution : 


( 5 - 3 ) 


. f fi; 1^1 f8 2^) and (5.3^) 

Steady state solutions of (5.1 ). (5.2 i aiu 

are obtainable by assuming 


Pq(o03 = " ^M-K+1 


(o03 = 0, 


solving successively the {M-K+2) 
Applying these conditions and solving 

' Ao 1. 

equations, the parameters 

|u^^ ^2 parameters are obtained as follows : 


p .(- l = -- 


1+ 


^ . ..4 1 


Ai 


il M 



Sterlllzalion is performed on attainment 


matrix is 


Case 11 : 

r^^rfnrmed on attainment 
i.e. sterilization is performea 

children. Parameter matrix Is t Xt> h 


of three surviving 



attainment of four surviving 


sterilization is made on 
matrix is I 


children 


7r - r I #. -t. ^0 j 

'■i"' H ^ 


, A^^Ai 

^3 

Xg Ai A-^ -] - I 

A/, /M, H, ■' 


TT. i tr 

M,H,_ ® 

^ - A^ A| A 2^ 

^ M, H^Hs " 

TT.. — /V Ai Av 


::: jV A| A%.A? A^t, j. 


^ 'I^ ■*■ TTm >«,>.3K ?.vAjA-M 


Case V: M-K+1 


e. sterilization is performed on attainment of six surviving 

4. • * f ^ # Ar A ^ A ^ A % A ^ V 

hlldren. Parameter matrix is ( ) 

^ r, A. .isMiia- 

L'*'7V 

.*. , * . %. ■ i 4 -.i ' ■ 


Xe X» AiAa 




X 9 ^ 
K ° 


'H' _ TT 

7VU r 

nr Af Xt- As . -jr- 

A# ^1 K A i..A*i.^ it„ 

lL,«* !!*■': '■ ■^"''" 1 '» j . £3 
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fii 


f^z 

A? ^2 Af 


Aj Av A3 Az Af A 


For presenting the effect of different levels of sterilisation 

at various parity levels of the couples 
sixth surviving children), we define the fertility 

mortality classes as follows : 


(from second to 


Fertility and Mortality 



steady state probability distributions of the number of surviving 
children classified by different levels of sterlization are 


given below 
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Table 

5.1 (F.M. 

Class I) 



Number of 
surviving 
children 

Level of 

2 

sterilization at parties 

3 4 5 

6 

0 

0.0500 

0.0154 

0.0056 

0.0025 

0.0014 

1 

0.2000 

0.0615 

0.0225 

0.0099 

0.0054 

2 

0.7500 

0.2308 

0.0845 

0.0373 

0.0203 

3 


0.6923 

0.2535 

0.1118 

0.0608 

4 

__ 


0.6338 

0.2795 

0.1520 

5 

_ 



0.5590 

0.3041 

6 





0.4561 


Table 5.2 (F.M. Class II) 


Number of 
surviving 
children 

2 

Level of 

3 

sterilization 

4 

at parities 

5 6 

0 

0.0323 

0.0071 

0.0017 

0.0005 

O.COOl 

1 

0,1613 

0.0355 

0.0086 

0.0023 

0.0007 

2 

0.8065 

0.1773 

0.0430 

0.0115 

0.0035 

3 


0.7801 

0.1893 

0.0508 

0,0152 

4 

- 


0.7573 

0.2032 

0.0608 

:5 " : 

- 

- 

- ■ 

0.7316 

0.2100 



Table 5.3 (F.M Class III) 


Level of sterilization at parities 


Number of 
surviving 


children 


Table 5.4 (F.M. 


Level of sterilization 


Number of 
surviving 
children 
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Table 5.5 (F.M. Glass V) 


Number of 
surviving 
children 

Level of sterilization at 
2 3 4 

parities 

5 

6 

0 

0.0187 

0.0028 

0.0005 

0.0001 

0.0000 

1 

0.1246 

0.0190 

0.0030 

0.0005 

0.0001 

2 

0.8567 

0.1304 

0.0207 

0.0034 

0.0006 

3 

- 

0.8478 

0.1346 

0.0223 

0.0038 

4 

• - 


0.8412 

0.1391 

0.0240 

5 




0.8346 

0.1439 

6 





0.8276 


Discussion ; 

An overall glance in the table 5.1 to 5.5 
reveals that the proportion of surviving children at all 
levels of fertility and mortality classes under consideration 
decreases consistently when sterilization levels are relaxed. 
For exarripie.Table I (F.M. Class I) shows that the proportion 
of couples having two surviving children in the long run 
is 75% when sterilization is made on attainment of two 
surviving children as against only 45% of the couples 
having six surviving children when sterilization level 
is relaxed to the extent of sterilizing couples only after 
having six surviving children. 
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The same trend is more or less revealed in all 
other tables, corresponding to different sets of fertility 
and mortality levels as may immediately be seen on inspection 


of 

the 

dia 

gonal 

elements of 

the 

matrices 

corresponding 

to 

each 

of 

the 

tables from I 

to 

V. However. 

a better 

insight 

into 

the 

variation of 

the 

proportion 

of 

surviving 


children as a result of changes in the level of parity 
under consideration for adopting a sterilization policy 
may be obtained when we consider the average no of surviving 
children under different F.M. set up. This is considered 
in the following table. 

Table 5.6 

Mean number of surviving children at different 
levels of sterilization 


F.M Class Level of 

. ■ 2 

sterilization 

3 

at number of surving 
4 5 

children 

6 

I 

1.7000 

2.6000 

3.4872 ■ 

4.3329 

5.0935 

II 

1.7743 

2.7304 

3.6917 

4.6485 

5.5957 

in 

1.8073 

2.7805 

3.7600 

4.7389 

5,7151 

IV 

1.8260 

2.8071 

3.7933 

4.7801 

5.7655 

V 

1.8380 

2.8232 

3.8130 

4.8036 

5.7935 


Table 5.6 provides, a clue for deciding the optimal 
level of sterilization depending on the number of surviving 
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children. For example, if the policy is to ensure at least 
one surving male child in the long run (as it is most 
customary in India) per couple then it may be seen that 
irrespective of the variations in the levels of the fertility 
and mortality the level of sterilization may optimally be 
made on the basis of three surviving children. Whereas 
if sterilization policy is relaxed to the extent of sterilizing 
couples only on the basis of more than three number of 
surviving children then one may expect a phenomenal 
increase in the growth rate of the population not consistent 
with the situation leading to adopt a sterilization policy 
motivated to cut down the growth rate of the population 
significantly. On the other hand if sterilization policy 

is made on the basis of only two surviving children then 
the condition of ensuring at least one male child in the 
long run per couple may not be realized although the 

same may fulfil the plan of reducing the growth rate 
in a fast increasing population. Such a policy may 

not be culturally and traditionally accepted by the people 
at large. 


CHAPTER-6 


ON A MARTINGALE APPROACH TO A PRORT EM 
ON STERILIZATION POLIC Y 

6.1 Introduction 

Suppose \ve have a problem of sterilization 

of mothers with number of surviving children, attained 

by addition of a birth at the starting point of the observation. 
These mothers may be considered homogeneous with respect 
to ages, i^ is a minimum number which makes the mother 
eligible for sterilization. 

However, if i^ is reduced to (i^ -1) first, 

then the mother concerned, is discarded for sterilization 

once for ail. Whereas those mothers whose number of surviving 

children move from i to (i +1) first (without passing 

o 0 

through any other state) is considered eligible for sterilization 
and they are sterilized at a random time on attainment 

of (i ^ +1) numbers of surviving children. Besides these 
two groups, the proportion of women who will remain in 

state number of surviving children throughout a period 

T from the starting point of the observation (the reference 

period TT which is considered large enough to cover all 

the three types of outcomes) will be discarded from the 

sterilization programme. 

The asymptotic estimator (for large T) of the 
above have been attempted to be obtained using Martingale 

stopping rule. 





6.2 Development of the Model 


."otations 


= Number of surviving children at any 

time t of the birth and death process 

with birth and death parameters : 

A i , A( i) and (i) respectively, 

i being the size of the population. 


for 


(ii) = P[X(t) = llx(O) = i^ ; X(j) = for 

>11 J< 0- TTo /Co ^ 

(iii) 0({T) = P[X(j} = i iX(0) = i^ for all j g(0,T)J. 


(iv) Pj^ (t) =Pp((t) =10 : k = 0,1,2, 


fie have 


P^(t,A ) =Pk(tKHA^ 




for k = 0,1,2,,....., 


subject to Pj^.j (t) = 0: for k=0. 


w? 


Let us define 


1 + ^1 + % 

. /fT ■ "^^2" 




4 

...18.2) 


Corresponding to the parameters ” Ad) and ^(i)i 


1= 1.2,3. 










Then it can immediately be shown that f(X(tM is a Martingale^ 
with respect to = I'lXCt) (a <r-field)j;^Janiri and Taylor ’ 

(1975)). Suppose we define a stopping time 

S_ - min.Ct ^ O.X(t) =i^,-l;X(t) = i^+TXf.O) = i^ ) 


c^given T ■ 


wne’^e 


T=Suj, 


,(6.3) 


represents the stopping time when X(t) 


1. or 


X(n = Iq +1 or X(t) = i Q 


for all (O.TJ. In this case 
= T, and constant time T is a Markov time which 

are considered as absorbing states of the process given 
X(o) = ip . Then by opti 
(Karlin and Taylor (1975)), 


X(o) = ip. Then by optional sampling theorem of Martingale 


EU(x(S^))). E(f(XO)) = E(f(ip) = t(ip (S-H 

since X(0) = ‘o ■ 

assumption. 

Then 

♦ Edd^)) = 






t6.r) 


f a,- 1) -{!-*> f6-B) 


1,1 h« qterlUzed and Cl-^) 

'hen r propcrtion will not be s e o 

o 


all be sterilized where 



the reference period in the sterilization programme 


Particular Case 


Therefore, the total proportion of mothers not to be sterilized 


where 


^ 1 

Ai A 2 

+ 

♦ 


A 1 A 2 

A. 

•+ 

^ A. A, 
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which 


Next, to obtain the expected stoppin 


time 


decides the expected time when policy decision is taken 

for sterilization, 100 (1- ) (1- ^ ) ^6 of mothers while 

1^1% mothers for sterilization programme. 
0 -* 


leaving 


We have 


Y(t) = g(X(t))-|^ ^ (X(u))(g(X(u)+l)-g(Xlu)) 


K(x(u)) (g(X(u))-g(X(u)-l)))du 


'f- (X(t)) for arbitrary function 
nf Y(tl exists. Further, if 


over 


ided the expectation 


g(i) be so chosen 


is a Martingale 


Y(t) = (g(X(t))-t) 


ECglX(t)3')-E(S ) 


E(Y(S )) = E(y{0)) 


by Martingale optional sampling theorem 
g(X(0)) = g(io ) 


Further 


E(g(X(S .p))) 


E(g(X(S )))-g(i ) 



we require the solution of the difference 


To obtain EfS ^ 
equation given by (6.11) 


Putting, g(i+l}-g(i) = ^ g(^) ~ vlij. 

We have the transformed equation of (6.11) as 


immediately get 


the solution 




and g(l) = 0 , g( 0 ) = 0 etc. 

In particular, for i = 3 . g(i^- 1 ) = I 


g(i q) = g( 3 ) 




" 4 - — l'f~h ^ 








X(S,j,))) = 


(i- + (i\ )f 

V ■ -f ('i^- I) - 

/ r -f ) - 0 ^ j) 


( r tn#) 

4 ( u - 0 -iiiLL--— -irt-- fi"<Wf'»+'J 
R‘,-')~4n-) — “^SCO 

+fi.) 

•■. .' .H':.. a ■. ■ „ i ■.. / i' ,V^' 




(6.7) respectively 


6.3 A Numerical lllustratioa 


05 and T » 3 year* 
















:: 


i 












The values of TT^ ^ oi , and E{S^) have been tabulated 


in table I, for 1-2, 3,4 to consider the suitability of the 
initial number i^ to be chosen for the sterilization programme. 


Table 6.1 


i 

0 

^0 

o< 

E(S.j.) 

{i-ir)(i-«) 

O 


(1) (2) 

(3) 

(4) 

(5) 

2 

0.0885 

0.7408 

1.476486 

0.2363 

3 

0.0868 

0.8106 

1.121893 

0.17296 

4 

0.0448 

0.9418 

0.154163 

0.0556 


The values of 7r„ 


0^ and E{S.j.) for different values of 


i, under hypothetical values of 

H H H ) 


3 ''* 4 ' ^ 


2' 3’ 

6.4 


Conclusion 


A glance in the Table 6.1 reveals the fact for 


i„ = 2. may be the suitablemost starting number of surviving 


children in this kind of sterilization programme under the 
hypothetical values of the parameters. This corresponds 
to about 24% of the mothers with 2 living children who 
may be sterilized. However, it will be interesting to employ 


may 



live 




live data providing estimates of X. and^j^, (i = 2,3,4) 
and this decide the policy of the sterilization programme. 












CHAPTER-7 

ON THE EVOLUTION OF DIRECT A ND INDIRECT STRATEGIES FOR 
THE REDUCTION OF CONSCEPTIQN RATE 
- A COUNTER MODEL APPROACH 


7.0. The efficacy of any family limitation strategy has 

to be judged primarily on its effect on the interval betv;een 
tv/o consecutive order of conceptions. Given that a conception 
has occured at the beginning of the interval, follov;ing 
which there exists an infecundable period (dead time in 
the parlance of Geiger Muller Counter Model theory): which 


is 

partially 

composed 

of a fixed 

gestation period (in case 

of 

a full 

term 

live 

birth) as 

well as a random period 

corresponding 

to 

the 

Post Partum 

Amenorrhoea (P.P.A). The 


period of Post Partum Amenorrhoea (P.P.A) is followed by 
the resumption of fecundable period which continues till 
a subsequent conception takes place. Hence statistically 
the interval between two consecutive order of conceptions 
can be conceived as the convolution of two probability 

distribution; the first one corresponds to the Infecundable 
period (a part of which is again fixed) and the other to 

the fecundable period during which a woman Is exposed 
to the risk of conceptions with the completion of the gestation 
period point chosen as origin. 





Development of the Model 


Assuming that the Gestation period is fixed and 

is of length Tr(= 9 months for a live birth) and the period 
of P.P.A. conforms to a negative exponential distribution 
with parameter A (Talwar (1964)) and the xvaiting time 

distribution for a subsequent conception follov-ing the resumption 
of the fecundable period is also a negative exponential 
with parameter (Singh (1964)) the waiting time distribution 
becomes a convolution of two negative exponential with parameter 
Aand fi respectively. The convoluted density function 


is given by 


where 


y; 0<X<0^ 

Al(t-y). 0 < /( 


Denoting the convolution of g 



Again shifting the origin to TT” to take account of the 


gestation prior to the period of Post Partum Amenorrhoea 


(which implies that the range of the random variable of 


the infecundable period plus the fecundable period 


lies between TT to oo i.e. TF <t <> oo), we have the expected 


time between two consecutive order of conceptions as 



Further assuming that the parameters A and /( vary 
different order of conceptions and denoting /\. and 


over 


the corresponding Poisson intensities (rate) for going from 


Infecundable period to Fecundable period and going from 


Fecundable period to Infecundable period respectively between 


order of conceptions 


) represents the interval between the (1-1) 


order of conceptions 


Methods for increasing the mean conception Interval V 
c.) i = 1,2,3,... can be realised In two different ways: 

(A) By reducing the intensity of conception {i=1.2.3 ) 

(direct Family planning (F.P) methods) (B) By reducing 
the intensity A of going from Infecundable period to fecundable 
period (1=1, 2, 3....). (Indirect method of Family Planning 
can be realised by Increasing the period of lactation etc.). 



is reduced to 


if the conception rate 


is considered to be 


where e 


a damping factor which may vary from one order of conception 
to other depending on Socio-Economic and cultural factors. 


ed interval under (A) is 


under (B) If the Intensity A; is reduced to 


J> 0 and,! - 1. v v 

then the increased interval under 





Finally under (A) and (B) implemented together the increased 


interval between (i-l)th and ith order of conceptions is 


given by 


Formulation of (7.7) is based on the rationale of 


interaction between direct and indirect 


methods of family limiation 


Numerical Exercise on the Basis of Data on South 


Indian Women 


The increase in the interval for order of conception 


in the light of the data based on a sample of women in 
South India (Srinivasan (1572)) details of which have been 


quoted by Mttra fl9B2) as follows 



of Indian Womens 


Interval between conceptions 


Average interval 
(in months) 


Variance of the 
intervaKin months) 


Monthly probability 
of conception! ) 


Non Susceptible 
Period (in months) 


stimates of n. (i = 1,2,3) are taken from the 
in the third row of the table 5.1. 


ctually the interval between two consecutive 
births which is the same as the waiting 
ween two consecutive order of conceptions 
aqsumDtlon of one to one correspondence 


Parameters 

First & 
Second 

Second & 
Third 

Third and 

Fourth 

1 

2 

3 

4 

Sample size 

297 

306 

245 





1,2.3) are obtainable by 


deducting the constant infecundable period due to gestation 


as 9 months and taking the reciprocal of the residual 


Further the assumptions relating to (A) and (B) 


v^ere made as follows 






Table 7.3 A 


Inter bir th Interval (1 -2) a s a result of^ 


different levels intervention 

by direct 

and 

indirect methods 

B I RTH I NTER V AL 1 N MO NTHS 

.BETWEEN 

1ST 

AND 2MD BIRTHS 


Degree of - 
Indirect 

Intervention 

Direct 

Without 

intervention 

With indirect 
intervention 

With 

direct 

interv- 

ention 

with direct 
& indirect 
intervention 

.95 

.90 

40.32 

40.65 

42.41 

42.73 

.95 

.05 

40.32 

40.65 

43.64 

43.96 

.95 

.80 

40.32 

40.65 

45.03 

45.35 

.95 

.75 

40.32 

40.65 

46.61 

46.93 

. 95 

.70 

40.32 

40,65 

48.41 

48.73 

.90 

.90 

40.32 

41.02 

42.41 

43.10 

.90 

.85 

40.32 

41.02 

43.64 

44.33 

. 90 

.80 

40.32 

41.02 

45.03 

45.72 

.90 

.75 

40.32 

41.02 

46.61 

47.29 

. 90 

.70 

40.32 

41.02 

48.41 

49.09 


Table 7.3 B 

showing the change in the Inter birth Interval (2-3) a s 
a result of different leve ls of intervention by direct 
and indirect method s 

BIRTH INTERVAL IM MONTHS BETWEEN 2nd and 3rd BiaiH 


Degree of 
Indirect 

Intervention 

Direct 

Without With Indirect 

Intervention intervention 

With 

direct 

Interve- 

ntion 

With direct 
and inter- 
vention 

0.95 

0.90 

42.83 

43 .23 

45.09 

45.49 

C.S5 

0.85 

42.83 

43.23 

46.43 

46.82 

0.95 

0.80 

42.83 

43.23 

47.93 

48.32 

0.95 

0.75 

42.83 

43.23 

49.64 

50.03 

0.95 

0.70 

42.83 

43.23 

51.60 

51.99 

0.90 

0.85 

42.83 

43.68 

45.80 

45.93^.^ 

0.90 

0.80 

42.83 

43.68 

46*09 

47 . 26 

0.90 

0.75 

42.83 

43.68 

47.43 

48* 76 

0.90 

0.70 

42.83 

43.68 

49.64 

50* 47 ' 

0.90 

0.70 

42.83 

43.68 

51.60 

52*43 



Table 7.3 C 


showing the change in the inter birth interval (3-4) as a 


result of different levels of intervention 
by Direct and Indirect Methods 


BIRTH INTERVAL IN MONTHS BETWEEN THIRD 4ND FOURTH BIRTHS 


Degree of 
Intervention 
Indirect Direct 


Without With Indirect Vlith direct 
Interve- Intervention Intervention 
ntion 


With direct & 

Indirect 

Intervention 


lance in the Tables 5.3A, 5-3B and 5.3C reveal 


An overall 


that inter conception (or inter birth) intervals substantially increases 


as a result of direct intervention whereas the indirect methods donot 


The same conclusion holds for all other 


appear to be much effective 


Further in almost all the entries the joint 


order of conceptions also 


linear effects of Direct and Indirect Interventions have been more or 


less found as the sum of the Direct and Indirect Interventions respectively 


howing the virtual absence of any Interaction between the two types 


if interventions 


GHAPTER-8 



ENERALISED P RQB ABI LI TY M ODBL_FOR M E ASUR J N 
IXTER GOMCEPTIVE DELAYS 


item of 


a community. A 


correct appraisal of Inter Conceptive 


delays should take 


into consideration (i) variation of 


fecundability 


different order of conceptions 


over 


(ii) variability of 


fecundability betv/een individuals at 


the same parity level 


The present exercise is devoted to develop system 


matically probabilit 


models concerning inter conceptive 


accounts of both the factors into 


Counter theory of type I with fixed dead 


consideration 


time has been employed to obtain the Palm Probabilities 


(viz. given that at a particular point of time a conception 
of a given order takes place the probability of the waiting 
time for the next conception after a fixed time t known 


obtaining the waiting t*roe for 


as Palm Probability of 


Estimation techniques concerning 


the next conception 


various parameters of the fertility process has also been 


attempted in the present exercise 



otation and Development of the Model 


The problem is to find out the waiting time of 
(i+lp*' conception given that the I conception has taken 
place at a fixed time, say 1=0. (Palm Probability); subject 
to the condition that following i*’’ conception at T=0. 

consisting of the gestation period 
follows, during which 


an infecundable period 
and the Post Partum Amenorrhoea 
the probability of a reconception is Zero 
the fixed dead time in a Counter Model type I) 


conception given 


The waiting time distribution of (if) 
that the i ^*^ conception has taken place 


has been 


derived in Chapter 


order 


rameters 


where 


of conceptions 


mnrp riaorously de derived 
result can more ngoiuuo y 


[lowever 


f. as follows 


ir ©vents (conceptions 

©vent ' { conception has 

y CtJ'. 


Probability 



Then the probability generating functions Ip-S-t-J 
,nd 0 (z.t) of V,^(t) andrfc(t) respectively are as follows 


Then, from Chapter ^ it follows 


We have 


V. (1980). Point 


and Isham 


(vide Cox 


London, p 


Chapman 


event given first one has happened 


time 


Then 

oa 

F^(t) PfT'h < 

Oa hi—, [- ^ 

Qfi-- — PT"' 




and Lf 




dt 


...(8.9) 


' . til" ' 

be the density function corresponding to the (n+1) arrival 
given that the first arrival is at T=0. Then 


Oo 


n 


FI 


-f ( t } ~ ^ y 

^ " I- z 

- 2,Chril f ^ \ 

cs I a -it ([-z)f 


k-t 2 


r o~') ! 

derived in Ghaptnr z. 


J= \ 


(8.10) 


Equating boih sides the coefficients of z ’.<e get f ^(t) 
viz. the waiting lime distribution of the 2nd everit given 
that the first has occurred at T=0. 

Hence 


f^{t) 


flk + l ){a+t ) 
(k+l) 


k+2 


(a+t ) 


k+2 


changing the parameters a.k to a ^ , k ^ for the (i+l) order 
of conception given that the i has taken place at T*0 


we have 



Further if there is a dead time iufecundable period 
after T=0 then the waiting time distribution is given by 






The mean & 2nd raw moment from the origin of the distribution 
(6.12) is as follows : 


-f- : — -t li- 


Hence the probability of (i+1) conception between (t-l) 
to t month given that i ’^'^ event has taken place at 

T=0 (1=1, 2, 3, . . • . ) subject to the condition that there is 

an infecundable period"^ follo^ving T=0. given by 




r(t) 


k,. 1 

a- (k.+l) 

1 1 ' 


(a.+t) 


:123:- 


:124:- 



Estimation of the Parameters 


Let us consider L = 1 i.e. the parameters of the 

waiting time distribution of the second conception g;iven 
that the first conception has taken place at T=0, subject 
to the condition that there is an infecundable period 


following T=0. 


Using the same data, of Srinivasan 


(1972) as referred by Mitra and Banerjee (1982) estimating 
equations for k and a^ are as follows ; 


345.8 + (36.1) = 1649.01 

k^(k^-l) k^ 


Direct solutions by the method of moments are given as 


2.93804 


a| = 56.11664. 


Similarly the using the mean and the second raw moment 
of the waiting time distribution between second and third 
order of conceptions as 38.9 and 1921.01 respectively 
and = 18.2. we have. 

= 2.90212 
a„ = 60.07392. 





ui^iiy ubLug /-i ^ = 38.6 andfl 2= 1870.06 months and 

^3" 6 months for the third and fourth order of conceptions 

e have 


3.01456 


60.29147 


The probability distributions of the women conceiving 
between (t-t+l) '^ months t = 0.1,2..... 59 for the (i-i+ll^^ 


order of conceptions (i 


1.2,3) assuming Lnfecundable 


1,2,3) following the 


conception is 


provided in table 6 respectively 



Table 8 


Probability dLstributioa of waiting time of conception 
by months between different consecutive order of 
conceptions measured from the date of expiry of the 
infecun_da_ble_j)eriod. 


Month 


Between First 
and second 
order of births 


Between Second 
and Third 
order of births 


Between Third 
and Fourth 
order of births 


0.06391 

0.05888 

0.05433 

0.05019 

0.04642 


0.06239 

0.05758 

0.05320 

0.04922 

0.04560 


0.06719 

0.06163 

0.05661 

0.05207 

0.04796 


0.04299 

0.03986 

0.03700 

0.03438 

0.03198 


0.04229 

0.03926 

0.03650 

0.03396 

0.03163 


0.04424 

0.04086 

0.03778 

0.03498 

0.03243 


0.02977 

0.02775 

0.02589 

0.02418 

0.02260 


0.02950 

0.02753 

0.02572 

0.02405 

0.02251 


0.03009 

0.02796 

0.02600 

0.02421 

0.02256 


0.02114 

0.01980 

0.01855 

0.01740 

0.01634 


0.02109 

0.01977 

0.01855 

0.01742 

0.01638 


0.02105 

0.01966 

0.01837 

0.01719 

0.01609 


0.01535 

0.01443 

0.01358 

0.01278 

0.01204 


0.01540 

0.01450 

0.01366 

0.01288 

0.01215 


0.01508 

0.01415 

0.01328 

0.01248 

0.01173 


0.01136 

0.01072 

0.01012 

0.00956 

0.00904 


0.01147 

0.01083 

0.01024 

0-00968 

0.00917 


0.01104 

0.01039 

0.00979 

0 . 0092 ^ 

0.00872 



0.00855 

0.00809 

0.00766 

0.00726 

0.00689 


0.00823 

0.00778 

0.00735 

0.00696 

0.00658 


0.00868 

0.00823 

0.00780 

0.00740 

0.00702 


0.00653 

0.00620 

0.00589 

0.00560 

0.00532 


0.00667 

0.00634 

0.00603 

0.00573 

0.00546 


0.00624 

0.00591 

0.00561 

0.00532 

0.00505 


0.00506 

0.00482 

0.00459 

0.00437 

0.00416 


0.00519 

0.00495 

0.00472 

0.00450 

0.00429 


0.00480 

0.00456 

0.00434 

0.00412 

0.00393 


0.00397 

0.00379 

0.00361 

0.00345 

0.00330 


0.00410 

0.00391 

0.00374 

0.00357 

0.00341 


0.00374 

0.00356 

0.00340 

0.00324 

0.00309 


0.00315 

0.00301 

0.00288 

0.00276 

0.00265 


0.00326 

0.00312 

0.00299 

0.00286 

0.00274 


0.00295 

0.00282 

0.00269 

0.00257 

0.00246 


0.00252 

0.00242 

0.00232 

0.00222 

0.00213 


0.00263 

0.00252 

0.00242 

0.00232 

0.00222 


0.00235 

0.00225 

0.00216 

0.00206 

0.00198 


0.06250 


0.06705 
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\ SU RVIVAL THE ORY APPROACH FOR OBTAINING THE DURATION 
OF _INTER _ LIVE BIRTH lOTERVALS' 


order of conception, (i=l, 2, 3 


vhere the origin of tine 


open 


of the fecundable period for the 


parity [cjafter the termination of the infecundable period 


(post partum amenorrhoea) corresponding to the (i-l) 


parity 


We assume a one to one correspondence between live birth 


and conception for the sake of simplicity 


Let the hazard rate following 


he termination 


n) be given as 


here ^ 0 and 
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If we denote the survival function R .(t) given by 


. o<±<ct. 

f » 


A |-£- -t,j £ ^ ct^ 

' ^ 

'A/t, 

■t t 

■fc. 


. £ . ■ ■ f . 

r Ctj'U) -Xd i*-- 


l.r -f > t, 


...( 9 .. 2 ) 


VarCT 


E(T.) = r R{t)dt 

I 0 ^ 

.) = /-f t R(t) dt - ( I R(t) dt)^ 


and the moments of 
origin is given by 

E(T^) 


k order (k = 1.2,3 


a(t)dt) 



Estimation of the Parameters 


be the ' mothers 


order of conception (and a live birth) in 


“•) respectively 


Also the times of having conception of the 1st, 2nd, ' 

^ J ^ bo --O / 

motherP^ recorded ' as ir > t • respectively 


We further denote the likelihood function of the sample as 


and L .respectively 


Then 


Precisely in a similar way 









on substitution of the estimate 


one can get the m.l.e. of 


Finally to estimate 


-:152: 


L 3 


= Q i- [- A; -tf - 


L -5 

"vxT" 


S: 

r -Hjt, -nj t (t;.-t,) - J 2 S. S.Ct,-tV-'''^s'') 


A,' ' 


- e 


/ 

I't^i ^'^ 32 .■* 


ii.- s/ 

_ M, t,- y>^(*2-V - 


O 


(9.10) 

on substitution of the estimated values of A^ and S. from 
(9.7) & (9.9) in (9.10) we get the r..Le. of On substitution 

I ' ■■ ' 

of the estimated values of A : .and A: . we get the estimates 

V ‘ ‘ th 

of the interlive birth interval between (i-1) parity and 

the i’^*’ prity (i = 1,2.3, ) given by 

00 

15~~BlTi) = / ^t) dt 

0 f. 

d.f 

t 


t. ^Ai' t- at 


0 




^^4- / ' G dt 

"i ® 


1 



,x; -Alt 


-A;i -a; £ ' ^ -* 




' '] ■^I"; 


>; i ; e 


. C/-t 


]■■ 


( 9 . 11 ) 


V^hile the variance 


of the waiting time is given by 




till the end of the. fertility 


By sufficiently extending t, 
span it is possible to make 
problem of censoring in the estimates procedure 


0 and thus avoid the 
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C/ 


^ari If j ~ 

A; 


If r, -Mt, 




■f t 




e“ ^ ^ ' (I t a; ye" ' { 








( 9 . 12 ) 


Vi’e begin v.’ith some hypothetical data, and using 

gs^ijn^Jc_sC 

the same to obtain the ' "m'arrfT^ i rue r v a 1 between 1st and 2nd 

A 


conceptio’i 

Let i 

h = 


A Numerical Illustration 



Then 


The normal confidence limit of the waiting time between 


1st to 2nd conception lies between 0.6054 to 4.0776 with 


more than 99 


level of 


Conclusion 


This simple method of maximum likelihood estimation 


over 


all parities and then expected waiting time over 


parities may be compared, to show the fertility 


behaviour 


The model admits 


generalisation in the case 


of the assumption of one to one correspondence between 


conception and live birth is dropped 


(iii) The model avoids the problem of censoring & complicated 


estimation techniques thereof 


(iv) The model admits further generalisation if individuali! 


variation is, fecundity within the same parity is 


considered; although the result remain: more complicated 



CHAPTER -10 


A GENERALISATION OF FREUND’S BIVARIATE EXPONENTIAL MODEL 
FOR INTERRELATION OF POST PARTUM’ AMENORRHOEA & LACTATION 


10.0 A wealth of earlier empirical studies have revealed 


that prolonged breast feeding following a child birth is 


helpful in postponing the date of resumption of fresh ovulatory 
cycle by prolonging the period of post partum amenorrhoea. 


Given that lactation has been discontinued at some time 


following a child birth at age t of the mother, the problem 


is to obtain the residual life time of the period of post 


problem 


partum amenorrhoea (or infecundable period) 


is of considerable importance in family limitation programme 


while studying how far the encouragement of the lactation 


process is helpful in preventing reconception. The problem is 
again similar to one in survival analysis of paired organs, 
viz., if out of two kidneys (parallel system), one fails 
the residual survival time of the second kidney is associated 
with higher hazard rate and hence its longevity is reduced. 


10.1 This kind of representation is provided in Freund's 
(1961) Bivariate exponential model given by 



Denoting lactation period (L.P.) by X and post partum 
amenorrhoea (P.P.A) by Y and £< $ being the intensities 


of expiring of L.P and P.P.A respectively subject to 


condition that on expiry of lactation the hazard rate for 

/ 

expiry of P.P.A. is changed from P to ^ and on expiry 


rat'' of expiry of lactation 


of P.P.A. first the hazard 


Usually p{ 


period is changed from 


The representation of (10.1) and (10.2) however do 
not take into account of the age of the mother. We replace 
o<, by c<lt) and by/3(t) where /3<(t) and /^(t) both 
being functions of ( ^ .t) and ( /3 .t) respectively without 

and being 


dependent on t 


^ -..{10 

oC^Cy 

li (t)X - i Bit) — 




A specific form say 



- ^ t 

H ^ 1 


and 

/ 3 ( t] -- t 


...(10.5) 


which correspond to Weibull failure rate. A specific 

example of the application (10.3) and (10.4) lies in finding 
out the conditional probability of expiry of L.P. and P.P.A. 
at the age of (t+x) and (t+y) respectively when (x 

< y); or P.P.A. and L.P. at the ages (t+y) and (t+x) 

respectively when y < x. 


10.2 Moments of the distribution 






e;X-t^y / 


t 

^ ^ n. I 




(t) (A(t)) 




J(^ ^ iiaW 

1 V 

(>i(t)i0Mf (xMiW)) 

(•<ft)t^('’^))^ (<[±\}\xltH(i(t\) 
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E(Xlt) 




Coefficient of ^ in the expansion of M(0j 


. © 2 I ' ^ 


flf'(t)+P(t) 

^'(t)(flr(t)+p(t) 


...( 10 . 8 ) 


Similarly 


E(X It) = 2 . _fe(t)(a^(t)4 8(t))fflt'(t)P(t)p(t) ^gL 

^(t)+p(t) (0f'(t))2 (eRtO^Ct)) 


...(10.9) 


E(Y|t) = 0 ((t)+B'(t) 

p' (t)(o((t)+^(t) 


...( 10 . 10 ) 


E(Y ^It) 


■■(■^(t)+^' (t) 


(6^(t)(^(t)+^(t) +1’ (t)^m"l-iD ? ). 

T|g'' Tt))2"Te^(t)+p(t) 


...( 10 . 11 ) 


E(XY|t) = Coefficient of 6 i^2 expansion of 

M(0^. 62 ! t) 


AWiLit) +ia)l ’it) l^i ^ ^ 

Wt)+p(^)^ eTttTFct) 


(10.12) 


Var(X|t) 


E(X^lt)-(E(Xlt)^ 


(Cii))^+ 2fl({t^lt)+ij6(t)i. 


(10.13) 


' 'r I . ; 'j .■ 






V ' -f'-; ' 
: • .-',5 \'..l ?' ‘ t- ; ■ ■ 
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10.3 Estimation of the 


e method of moments 


of the mother at the child birth 


Given 


E(yjt), Var(Xjt) and VarfYjt 


at different levels of t 


One may combine the estimates assuming that the 


t of the mother 


Then using (10.15), one can get the correlation 


XY /Var(Yi t) 
/Var(Xjt) 




=ind 


Y = E{Y|t) + ^XY rx-EfXitl 

-/Var(x| t) 

^^yyV arTxft) 


(10.16) 


X = E(X t) + 


\/var(y| 


(Y-E(Y t) 


...(10.17) 


t) 


(10.16) will thus predict the period of P.P.A. given L.P. 
when L.P. < P.P.A. 

(10.17) will predict the period of L.P. given P.P.A. when 
P.P.A. < L.P. 


A numerical illustratio n of the technique 


The purpose of the Freund's model is to predict the 
Post Partum amenorrhoea (PPA) given the lactation period 
and lactation period given the PPA. These may be donn 
by using the results (10.16) and (10.17) respectively by 


the estimates of the. parameters tc. #((in flf(t) 

«C-1 - 


t ), $ Cn B (t) = p t' P respectively 


0 ^ -j^,k.$onXlglAK/ . I . 

by '.the method of moments given the sample mean E(Xjt)), 

EfYl’t), Varfxlt) and Var(Ylt) at different levels of t. 


Let us take the hypothetical example giving 


E(X|t) = 1.89 


2.60 


A 


Var(X|t) = 3.08 


A 

E(Ylt) 

A 

Var(Y|t) = 5.65 at any t 
















using the equations (10.8), (10.10) 

the estimates of the parameters ^ ^ 

by successive iterations as follows 


are obtained 


Putting these estimate 


Cor(X,Yit) ’ = 15) 


on substitution 


Var (X;t' ^ Var(Y|t) and 

Gorr('x.Ylt) in (10.16) and (10.17) 


given 


respectively we get the prediction 


Y = 2.89 months 


Conclusion 


the estimated prediction 


Precisely 


However 


Xdactation period) given 


"■ is ' r'" 

'A. , S' 4 i ■ : » - 

: f ' ' 'A 
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examples are'^^ypothetical to illustrate further application 


of the model based on real data. It Is in fact, difficult 
to obtain the empirical joint distribution of P.P.A. and 


lactation period.. Besides that the problem of censorinp; _ 

ph^lfdb 

also there even a follow up ^aTtev a long gap.^' KJnder 


is 


these difficulties it is possible to employ the model to predict 
J 


P.P.A. given lactation and vice and versa. However, the 
real efficacy of the methodology, can be understood when 
more extensive live data are obtained. 


>l‘r' '■ 





CHAPTER-11 


CONCLUSIONS 


The studies in this dissertation pertain to stochastic 


modelling of population in the area of fertility analysis 


and family planning programmes. The study begins with 


a technique based on Palm Probability for the comparison 


of fertility statuses of tv;o otherwise socio-economically identical 


the age of effective 


The idea is to examine the policy of reduction of fertility 


by postponement of marriage. The solution which has been 


derived is purely methodological one applied on hypothetical 


Nevertheless, it is believed that the methodology 


data 


developed may be tested on a more exhaustive set of real 


data to bring out more meaningful result. Similarly in 


chapter four a model based on Palm Probability, so developed 


has been employed to obtain the waiting time distribution 


aiting time distribution 


of first conception from marriage and 


between consecutive order of conceptions 


The traditional model building exercise has been 


Palm Probabilistic technique 


modified here by using 


But the application of the technique has been done while 


keeping restriction of infinite marital exposure. In this 

connection comparison of the results based upon the approach 


of infinite marital exposure with that of finite marital 



exposure (Sresthn and Biswas (1985)) reveals that the gain 
in precision in the estimate of interconception intervals 
for lower order of conceptions is almost negligible. However, 
for higher order conceptions further investigations are necessary 


Attemps have been made in chapter five to develop 
a Multistage Markov Chain Model based on a density dependent 
birth and death process to obtain the asymtotic probability 
of having fixed number of surviving children based upon 
the same distribution. Attempts have also been made to 
decide an optimal sterilization policy depending on the 


number of surviving children 


One major indication in this study, although a methodol- 
ogical one based on hypothetically chosen fertility and 
mortality rates, is that if the policy is to ensure at least 
one surving male child in the long run (i.e. it is essentially 
desirable lower limit in India) then the result shows that 

irrespective of the variation in the levels of fertility and 

mortality the level of sterilization may be made on the 

basis of three surviving children. If policy is relaxed 

to the extent of sterilizing mothers only on the basis 
of more than three surviving children, then one may expect 
higher growth rate which may not be consistent With 

the policy of reduction of growth rate. Finally, if ’ an 
intermediate action based on the policy of two surviving 

children is adopted then the condition of ensuring one surviving 




male child in the long run per couple may only be realised 



with very low probability 


Although the same may fulfil 


the plan of reducing the growth rate, it may not be culturally 


accepted by the people at large 


In chapter six several martingales have been constructed 


on the continuous parameter birth and death process to 


obtain the stopping time for sterilizing a mother and also 


estimating the proportion of eligible mothers for sterilization 


in future, given that each of them has a fixed number 


of children at present 


This study is based on hypothetical data. Nevertheless 


this kind of methodological investigation is useful in answering 


two basic problems 


one 


sterilization programme. The percentage of mothers obtained 


for the sterilization programme based on the choice of different 


number of surviving children is indicative of the increase 


or decrease of the probable cost if we switch over a more 


liberal or more conservative type of policy while adopting 


The second one relates to 


the sterilization programme 


one can indirectly calculate the 


time i.e 


of the mother at the time of sterilization connected 


with various sterilization programme. Further investigation 
of these two items of consideration may be of crucial importance 


The appliction of Martingale 


in our national planning 


theory in such problems with live data would indeed be 
interesting. 



Chapter seven deals with the evolution of direct 


and indirect strategies for reduction of conception rate 


based on an approach of Geiger Muller Counter Model 


The direct strategies correspond to that of mechanical 


methods of family planning measures, whereas the indirect 


trategies lead to lengthening the period of lactational 


amenorrhoea and abstinence period (i.e. temporary 


seperation of husbands and wives from conjugal life 


following a child birth) as it is necessitated by social 


The methodology has 


customs in several communities 


been applied on Srinivasan's data (1972) based on 


the experience of a sample of South Indian Women but 


the rates of conception as it may exist under the assumption 


direct and indirect strategies of family limitation 


programmes on both are again based on hypothetical 


data. As such the investigation again becomes a methodolo- 
gical one. Nevertheless such investigation may highlight 
the rate of both direct and indirect strategies of family 


The relative 


limitation in reducing the birth rate 


efficacies of the factors are to be taken into consideration 


for evolving suitable family planning programme. 

In Chapter eight a generalised probability 
model for measuring inter-conceptive delays have been 
attempted. Inter-conceptive delays inthis study haw 
been taken as the function of order of conception as 


well as individual fecundability affected by biological 





and social condition. The non-susceptible period 
between conceptions have been accounted bv the dead 


time of the counter model and the analysis has been 


attempted to be more precise by considering the probability 
of a conception based upon previous fertility behaviour. 
Further application with real data may be fruitful 


for this kind of methodological exercise 


Chapter nine deals with the survival theory 


approach for obtaining the inter live birth interval 


based upon a survival function which has been chosen 


to be quite flexible i.e 


it can take into account 


type of fertility behaviour, viz 


increasing fertility 


decreasing fertility with age or a fertility 


function showing increasing as well as decreasing fertility 


behaviour at different age segments 


The estimation of the 


of the model 


by the method of maximum likelihood is straight forward 


This has been illustrated while taking the hypothetical 


data concerning the first & the second order of conceptions 


Hov;ever the same can be easily generalised for higher 


order of conceptions 


been considered in this kind of approach is that it 


avoids the problem of censoring. Even one to one 


correspondence between live birth and conception as 


essential for building up a generalised model. However 



with such realistic assumptions mathematical complications 


would be quite considerable but this would be indeed 


a worth taking exercise. Further exploration of the 


exact status of the maximum likelihood estimator in 


this case may also be very useful 


Finally chapter ten deals v;ith a study of the 
role of lactation period in affecting the period of Post-partum 
amenorrhoea. The novelty of the methodological approach 


lies in generalizing Freund's model which is a bivariate 


exponential model with certain degree of dependence 


between the time of failure of two components or two 


sources of failure subject to the age dependence of 


the failures: where the failures are defined as 


period of completion of the Post Partum Amenorrhoea 


and stoppage of lactation period of mothers respectively 


Post partum amenorrhoea is a natural factor 


but with the cessation of the same the lactation period 


discontinued with greater intensity interestingly 


comparable with that of the time of failure of a second 


component given that the first component has already 
failed in a bilateral parallel system. Further^ demographical ly 
speaking the problem has become more complicated because 
of the dependence of the two factors being influenced 
by the age of the mother. A generalisation of Freund’s 
model has been found to be quite fruitful tool to analyse 



However the limitation of the 


hypothetical data 


A meaningful conclusion on the 


applicability of Freund's model as well as the result 


based on the same would be an useful item of information 


based upon exhaustively collected real data 


The application of Palm Probability, Counter 


Model 


Rene%val Theory and other stochastic process 


oriented techniques as Martingale Theory are some new 


methodological 


applications for arriving at solution 


of ferility analysis 


However the limitation 


ofthe 


study is that all the studies are methodological and 


all solutions have been evolved only on theoretical 


plane but the applications of the nev/ly evolved techniques 


and their justification will much depend naturally on 


more exhaustive type of real data. This is the motivation 


of making the methodological exercise for a fruitful 


application in the area of fertility investigation 




APPENDIX 


To show 


A' (T-tJ (k+1) (k+2)( 


(k+3) 


given that 


(k-»n) 



(k+n) 


a+A(T-t , ) 


a+A {T- t, ) 


(k+n) „k+l 


n^(k+l) (k+2) 





(k+l)(A{T-t,) 


(k+l)(k+2){ 


a+A(T-t- ) 


(l+(k+2)( 


a+A(T-t , ) 


+ (k+2)( 


a+A{T-t, ) 


a +A (T—t ^ ) 


(k+2) 


a+A{T-t,) 


(k+3) 



(k+1) A'(T-t,) A (T-t.) 


(a+A(T-t,)) 


A(T-t.) „ 

A- (T-tj) (kn)(k.2) 


A(T-tj) ■'t^+2) 
a+A(T-t, 


A(T-tj) 

l^^a+Ad-t, j 


Ck+l)a 



Papers prepared from the contents of this thesis 


On a Palm Probabilistic Approach for the Comparison 


of Cohort Fertility; Demography India, Vol. II. No. 1, 1982 


On Probability for Obtaining Inter-Arrival Time Distribution 


in Weighted Poisson Process; Calcutta Statistical Association 


Bulletin, Vol. 32. Nos. 125-126, March and June, 1983 


On the application of Palm Probability for Obtaining 


Waiting Time Distribution Between the First and Higher 


Order Conceptions; Proceedings of the XIII International 


Conference on Stochastic Processes and their Applications 


Demography India vo. 18, No. 2 page 277-312 


A Multistate Markov Chain Model for Evaluating ,a Sterilizationn 


Policy, Biometrical Journal vol. 29, No. 1, Page 57-64 


On a Martingale approach to a problem of Sterilization 


policy 


Stochastic Modelling and Decision Making 


University of Delhi, March 1986 


Khanna Publishers 


New Delhi 




BIBLIOGRAPHY 


1. Basu, D. (1955) : A note on the structure of a stochastic 
model considered by V.M.Dandekar ; Sankhya Series B. 
15. 251-252. 

2. Bhat U.IJ. (1971) : Elements of Applied Stochastic Processes 

- John Wiley & Sons. New York. 

3. Bisv^as, S. (1963) : A study of amenorrhoea after child 

birth and its relationship to lactation period; 

Indian Journal of Public Health. 7. 9-14. 

4. Biswas. S.(1971) : On a probability model of the 

waiting time of conception based on Censored Sampling: 
Proceedings of the All India Seminar on Demography 
and Statistics, Banaras Hindu University, page 257-264. 

5. Bisv/as, S. (1973) : A note on the generalisation of 

William Brass Model : Demography 1C, 459-467. 

6. Biswas, S. (1973) : Abstinence, Post Partum Amenorrhoea 

and interpregnancy interval - Demography India, Vo.2, 
page 203-211. 

7. Biswas, S. (1975) : On a more generalised probability 

model of the waiting time of conception based on censored 
sampling from a mixed population; Sankhya Series B, 
37, 343-354. 

8. Biswas, S. (1981) : On the extension of some results 

of counter models with Poisson inputs and their 

applications ; J. Indian Statist. Assoc. 18. 45-53. 

9. ^^^^ B^^^^ S. and Pachal T.k. (1982): On a Palm Probabilistic 

approach for the comparison of cohort fertility - Demography 
XI, No. 1, page 135-140. 

10. Biswas, S. and Pachal T.K. (1983) : On the application 

of Palm Probability for obtaining inter-arrival time 

distribution in weighted Poisson Process ; Cal. Statist. 
Assoc. Bulletin, 32, 111-115. 

11. Biswas, S. and Pachal T.K. (1984) : On the application of 

Palm Probability for obtaining the waiting time distribution 
between the first and the higher order of conceptiohs 

- Demography India, vol. 8, No. 2, page 277-312, 






David and Moeschberger M.L. (1978) 
of Competing Risks Monograph No. 39, 
Griffin and Co. Ltd., London. 


The Theory 
Charles 


Dharmadhikari, S.^, (1964) : A generalization of a 

stochastic model considered by V.M. Dandekar; Sankhya 
Series A, 26, 31-38. 


Biswas. S. and Pachal, T.K. (1986) : A multistate 
Markov Chain model for evaluating a Sterilization 
policy - Biometrical Journal, voL 29, page 57-69.^^^^ ^ ^^ ^^^^^^ ^^^^ 

Biswas, S. and Pachal. T.K. (1990) : On a Martingale 

approach to a Sterilization Policy - Proceedings of the 
U.G.C. sponsored Seminar on 'Stochastic Modelling 


and Decision Making', University of Delhi. 
Khama Publishers, New Delhi, page 227-235. 


March 1986, 


Bogue (1980) 
Community and 
Dhicago. 


Comparative birth interval analysis: 
Family Studies Center, University of 


Bongaarts and Potter (1983) 
Behaviour: Academic Press. 


Fertility, Biology and 


Brass, W. (1958) : The distribution of births in human 

population: Population Studies 12, 51-72. 

Braun, H.I. (1980) : Regression-like analysis of birth 

interval sequences; Demography 17, 207-223- 


Cox. D.R. and Isham, V. (1980) 
Chapman Hall, London. 


Point Processes: 


Chiang, C.L. ( 1968) : Introduction to Stochastic 

Processes in Biostabilities - John Wiley and Sons. 

Chiang, C.L. (1979) : An introduction to Stochastic 

processes and their application - Robert E, Krieyer 
Publishing Company, Hanlington, New York. 

Dandekar, V.M. (1955) : Certain modified forms of 
Binomial and Poisson Distributions; Sankhya Series 
b, 15, 237-250. 

Das Gupta Prithwis (1972) : On two-sex models leading 

to stable populations; Theoretical Population Biology 3, 
358-375. 



D'Souza. S. (1974) : Closed birth intervals 

analytic study; Sterling Press. 


A data 


Deund J.E.(1961) : A bivariate extension of the exponential 
distribution - J.A.S.A. vol. 56, page 971-77. 


Gini (1924) : Quoted in losifescu M. and Tautu P. 

(1973) Stochastic Processes and applications in biology 
and medicine, II Models, Biomathematics 3; Springer 
- Verlag, New York- 


Gross, A.J. and Clark A.V. (1975) .- Survival Distribution 
- Reliability application in Biomedical Sciences - 
John Wiley, New York. 


Henry, L. (1953) : Fundaments theoriques des measures 

de la fecondite naturelle; Revue de I'lnstitut International 
de Statistique 21, 135-151. 


Henry, L. (1957) : Fecondite ef famille. Modeles 

mathematiques (I); Population, 12, 413-444. 


Henry, L. (1961a) : Fecondite ef famille. Modeles 
mathematiques (II): Population 16, 27-28. and Population 
16, 261-282. 


Henry, L. (1961b) : La fecondite naturalle observation 

theoric, resultats; Population 16. 625-636. 


Henry, L. (1961c) ; Some data on natural fertility 

Eugenics Quarterly 8, 81-91. 


losifescu, M. and Tautu, P. (1973) : Stochastic processes 
and applications in biology and medicine I and II 
Models : Biomathematics 3; Springer - Verlag, New 
York. 


Karlin, S. and Taylor, H.U. (1975) 
in Stochastic Pro(^sses, Second Edition 
New York. 


A First Course 
Academic Press, 


Keyfitz, N. (1966) : Sampling variance of Demographic 

Characteristics - Human Biology, vol. 18, page 22, 


Keyfitz. N. (1967) ; Reconciliation of Population Models" 

Matrix Integral Equation and Partial fraction - Journal 
of Royal Statistical Stociety, Series A, vol, .136, 
page 61. ' V' 

Keyfitz, N. (1976) : Introduction , to the Mathematics 

of population with revisions Addlson-Wtsley Publishing 
Company - Reading, .Massachusetts,' lit.. 



Lawless J.F. (1982) 
for life time data. J 


Statistical Models and Methods 
y, New York. 


Khintchine, A.T. (1960) : Mathematical methods in 

the theory of queueing (Translated by D.M. Andrews 
and M.R. Quinouille) ; Charles Griffin, London. 


Menken, J. and Sheps, M. (1977) ; Mathematical models 

of conception and birth. 


Mitra, S. and Banerjee, S.N. (1982) : An analysis 

of intervals between successive orders of birth; Demography 
India II, 115-123. 


Neyman, J. (1949) ; On the problem of estimating 
the number of schools of fish; University of California, 
Publication in Statistics 1, 21-36. 


Pathak, K.B. (1966) : A probability distribution for 

the number of conceptions : Sankhya, Series B, 28, 
213-219. 


Pathak, K.B. (1970) ; A time dependant distribution 

for the number of conceptions: Artha Vijanana, 12, 
429-435. 


Perrin, E.B. and Sheps, M.C. (1964) : Human reproduction 
A stochastic process; Biometrics 20, 28-45. 


Potter, R.G., Sagi, P.C. and Westoff, C.F. (1962) 
Improvement of contraception during the course of 
marriage; Population Studies 16, 160-174. 


Potter, R.G. (1970) : Births averted by contraception 

an approach through renewal theory: Theoretics 

Population Biology, 4, 251-272, 


: A renewal process, relation to type 
counter models - Annals of Mathematical 
, pp. 737-754. 


Pyke R. (1969) 

I and type II 
Statistics vol. 29 


Fertility differentials 


India 


Registrar General 
in India, 1972. 


G. and Hobcraft, J.N. (1980) : Illustrative 

Life table analysis of birth intervals In 
World Fertility Survey Report 16. ; 


Rodriguez 

analysis 

Colombia 



Ross, J.A. 
model for 
35, 439-454 


and 

birth 


oneps, M. C. aind Menken, 
of birth intervals accordir 
Theoretical Population Biology 

Sheps, M.C. and Perrin, E. 
of birth intervals under a 
Models; Population StnrfieQ 17 


J-A. C1972) 
f to the 
3, 1-26. 


: Distribution 
sampling frame; 


bJieps, M^G. end Menken 
studying birth rates g 
in reproductive parameters; 

Singh, S.N. ( 1963) : r 
in the number of births 


J.A. (1971) : A 

given time dependent 
j; Biometrics 27, 325-343 

Probability models for the 
J per couple; J.A. S.A, 58, ; 


variation 

721-727. 


Singh, S.N. (1964a) : On th 

Sankhya Series B 25, Part-1 and 2 

Singh, S.N. (1964b) : A proba 

fertility; Sankhya Series B, 26, 89 


couple 


bingn, b.N. and Pathak, K.B. (1968) 
of the number of conception; 
Journal, Ban aras Hindu University, 1, 

Singh, S.N. and Bhattacharya, B.N. 1 
probability distribution for couple 
26, 33-39. 


and Bhadurl, T, (1971) : On the 

rtum amenorrhea; Proceedings of All 
Demography and Statistics. Banaras 
Varanasi, India. 


; Counter 

Camb.,PhiL 









''M 




■ 


Smith, W.L. (1958) : Renewal Theory and its Ramifications 
- J.R.S.S. Series B, vol. 20, pp. 243-302. 


Srinivasan, K. (1968) : A set of analytical models 

for the study of open birth intervals; Demography 
5, 34-44. 


Srinivasan, K. (1972) : Analytical model for two types 

of birth intervals with applications to Indian Population; 
St. Joseph's Press, Trivandrum-695014. 


Takacks, L. (1960) : Stochastic Process (Problems and 

solutions) ; John Wiley, New Delhi. 


Talwar, P.P. (1965) : A bio-statistical model on 

the mechanism underlying the spacing between two 
consecutive births and its applications to the study 
of the effect of contraception: Paper No. 2 in the Seminar 
on new approaches to the use of mathematical models 
in demographic analysis at the Demographic Training 
and Research Centre, Bombay. 


Yadava, R.G. (1966) : A prdbability model for the 

number of births; Seminar volume in Statistics, Banaras 
Hindu University, Varanasi, India. 


■Additional Referenceat 




1, Sehgal V.K, Pachal T.K, ( 1994)1 On a generalised age dependent aaodel 
for the estimation of mean Survival tise in* the presence of tijo oompeting- 
risks- Under publication in International Journal of Systems Soieoce* 

2, Biswas S and Sehgal V.K. ( 1982)1 On a generalised proUaMlity model 
for estimating the proportion of female popilatiott at diff«P«Btt levels 

of fecundity based on censored sampling from a mixed paimlat|,Ottp- InterMational| 

journal of Systems Science, Vol.l8»No«I0 »Fage Z909*«19J7*k r \ ; 

3, Srestha Oanga and Biswas 3 ( 1985)1 International report tO/8^/$0‘IGfP 

Trieste,! taly.. "'-i / ,■ 








